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ABSTRACT 


This dissertation consists of a theoretical and experimental study of the 
performance limits of a free-space direct detection optical communication 
system using a semiconductor laser diode as the optical transmitter and a 
silicon avalanche photodiode (APD) as the receiver photodetector. Optical 
systems using these components are under • extensive consideration as 
replacements for microwave satellite communication links. Due to the peak 
and average optical power output constraints inherent in semiconductor 
laser diodes, optical pulse position modulation (PPM) was chosen as the 
most appropriate signal format. 

An experimental system was constructed that used an AlGaAs 
(Hitachi-HL8314E, X=833nm) semiconductor laser diode as the transmitter 
and a silicon avalanche photodiode (RCA-C30902S) photodetector. The sys- 
tem used Q=4 PPM signaling at a source data rate of 25 megabits per 
second. The PPM signal format requires regeneration of PPM slot clock 
and word clock waveforms in the receiver. A nearly exact computational 
procedure was developed to compute receiver bit error rate without using 
the Gaussian approximation. A transition detector slot clock recovery sys- 
tem using a phase lock loop was developed and implemented. A novel word 
clock recovery system was also developed in which a phase lock loop was 
used to track the random occurrence of back-to-back PPM pulse pairs on 
the received data stream. This enabled the receiver to establish PPM word 
boundary synchronization without the presence of specially inserted syn- 
chronization patterns in the transmitted data. 



It was found that the results of the nearly exact computational pro- 
cedure agreed well with actual measurements of receiver performance. Use 
of the widely accepted Gaussian model for the APD output photocurrent 
was shown to underestimate the value of APD average gain that minimized 
the bit error rate of a maximum likelihood receiver structure and to sub- 
stantially overestimate the actual bit error rate of the receiver when the 
average number of detected background noise photons was less than one per 
PPM slot time. The maximum likelihood receiver achieved a bit error rate 
of 10 -6 at received optical signal levels which corresponded to 55 average 
detected signal photons per bit. This is best receiver sensitivity yet reported 
for a direct detection system of this type. The transition detector slot clock 
recovery system was simple to implement and yet caused no measurable 
penalty in receiver sensitivity when the phase error was less than one per- 
cent of a slot clock cycle. This occurred when the noise bandwidth of the 
phase lock loop was less than 2X1CT 4 of the slot clock frequency. The new 
word clock recovery scheme was shown to perform well in initial synchroni- 
zation. The entire PPM slot and word timing recovery subsystem could 
acquire and maintain receiver synchronization at received optical signal lev- 
els corresponding to more than 15 detected signal photons per transmitted 
information bit. 
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Chapter 1. Introduction 

Fiber optic guide laser communication systems have revolutionized ter- 
restrial long distance data transmission. Advances in the technologies of 
semiconductor laser diodes and silicon avalanche photodiode (APD) photo- 
detectors have now made it possible to realistically consider replacing exist- 
ing microwave satellite communication systems with optical communication 
systems. Semiconductor laser diodes are small and highly efficient devices 
which require only low voltage and low power (few watts) electrical power 
sources to generate tens to hundreds of milliwatts of optical output power. 
Because of the much shorter wavelength of the optical frequency radiation 
(1 micron versus 1 millimeter), very large antenna gains (~100dB) can be 
obtained with relatively small sized (10 cm diameter) optical components. 
As a result, high data rate (100-500 megabits/second) free space optical 
communication system have become highly competitive with comparable 
microwave systems [l, 2, 3, 4j. 

The main objectives of this dissertation were to conduct a theoretical 
as well as experimental investigation of the performance limits and receiver 
timing recovery of a direct detection optical communication system which 
used currently available semiconductor laser transmitter and silicon APD 
photodetectors. The net results of this work was a state-of-the-art system 
that used Q=4 pulse position modulation (PPM) signaling at a source data 
rate of 25 megabits per second. The system achieved a receiver bit error 
rate (BER) of 10 -6 at received optical signal levels which corresponded to an 
average of 55 detected signal photons per source date bit. This is within 10 
dB of the absolute (quantum) limit to system performance and is the best 
performance yet reported for this type of system. 
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A theoretical analysis of the system performance revealed that stan- 
dard Gaussian models for the APD output photoelectrons were insufficient 
for proper system design. A more accurate model was used which required 
the development of new numerical computational procedures efficient 
enough to allow receiver BER to be computed without the use of excessive 
amounts of computer time. Novel PPM slot and word boundary timing 
recovery subsystems were also developed. Their performance were 

mathematically characterized and shown to agree well with experimental 
measurements. Use of the recovered clocks caused no measurable penalty in 
receiver sensitivity when compared to a receiver with common 
transmitter/receiver clocks. The net result was the development and 
laboratory demonstration of a completely self synchronized PPM receiver 
that could acquire and maintain timing synchronization with the received 
PPM data stream at received optical pulse energies as low as 15 average 
detected signal photons per source date bit (30 photons per light pulse). 
The remander of this introductory chapter further describes the receiver 
design and analysis and gives explanations for the various design choices 
made. 

A direct detection optical communication system consists of a photo- 
detector which directly converts the received optical signal intensity into a 
photocurrent. The output of the photodetector can be modeled as a Poisson 
counting process in which the average counting rate is proportional to the 
received optical field intensity [5]. The optical signal can only be intensity 
modulated in a direct detection system. 

Another type of optical communication system is called a coherent 
detection system. The receiver in a coherent detection system must contain 
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a local oscillator which is phase coherent with the received signal optical 
field. The local oscillator laser is superposed with the signal optical field 
and a photodetector is used as the nonlinear device (square law) to generate 
an intermediate frequency (IF) component. The IF component is propor- 
tional to the input signal field strength rather than intensity, and therefore, 
the optical signal can be either amplitude, frequency, or phase modulated. 
The performance comparisons between coherent detection and direct detec- 
tion systems are given in [6, 7], In general, direct detection systems are far 
simpler to implement. 

Direct detection optical communication systems are sensitive to back- 
ground radiation because the receiver noise bandwidth is determined by the 
optical filter preceding the photodetector which is typically on the order of 
a few angstroms, or hundreds of gigahertz. Coherent detection systems, on 
the other hand, are relatively immune to background radiation because the 
receiver noise bandwidth is determined by the IF filter bandwidth which is 
usually close to the source data rate. The performance of direct detection 
systems is also limited by amplifier thermal noise. A coherent detection sys- 
tem, on the other hand, can override amplifier thermal noise by increasing 
the local oscillator laser power, since the IF signal component is propor- 
tional to the amplitude of the local oscillator. The performance of coherent 
detection systems are limited by the amplitude and phase noise of the local 
oscillator laser. A coherent detection system requires a very complicated 
optical phase tracking loop in order to generate a coherent local oscillator. 
The local oscillator laser field and the signal optical field have to be aligned 
perfectly in spatial mode and polarization in order to generate the IF com- 
ponent properly. Although the technologies of coherent optical communica- 
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tion are advancing rapidly, semiconductor lasers available today still have 
excessively high levels of phase and frequency noise. The stability of laser 
devices, the complexity of the receiver, and the reliability of the entire sys- 
tem are still problems to be overcome in order to fully exploit the potential 
of coherent communication systems [8j. Direct detection optical communi- 
cation systems, on the other hand, are much simpler, less expensive and 
more mature. They have been widely used in terrestrial fiber optic commun- 
ication systems and will surely play an important role in free space optical 
communications. 

A PPM signal format can be used to improve performance of a direct 
detection optical communication system when the optical transmitter is 
average power limited, rather than peak power limited [9]. In PPM signal- 
ing, a group of L binary source bits are transmitted as a PPM word which 
consists of a single light pulse positioned in one of the Q=2 L time slots, as 
shown in Figure 1.1. If the binary source data rate is 1/T B , with T B the 
pulse width of the corresponding input binary on-off-keyed (OOK) signal, 
the PPM pulse width, T s , is given by QT S = LT B , or T s == (log 2 Q/Q)T B . 
The major advantage of PPM signaling format is that it concentrates the 
limited average transmitter laser power into fewer narrower, but higher 
peak power, PPM pulses. The system performance, as measured by the 
probability of correct detection, improves because the signal becomes more 
distinguishable from the background and circuit noise. Since the signal 
bandwidth of the system is proportional to the reciprocal of the pulse 
width, PPM signaling requires Q/log 2 Q times larger bandwidth as that of 
OOK signaling. However, unlike optic fiber systems, the bandwidth of a 
free space optical communication system is limited only by the electronics 


Binary 

bit 1 bit 2 



bit L 


Source Data t J 




PPM Signal 
(Q=2 L ) 




QT,=LT 


Figure 1.1. PPM signal format 
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rather than by the light dispersion properties of the transmission media. 
AlGaAs laser diodes as optical transmitters are average power limited for 
Q<8 PPM signaling (10, 11]. Therefore, low order PPM signaling is suit- 
able for free space direct detection optical communication systems when 
AlGaAs laser diode transmitters are used. 

The ultimate performance of an optical communication system is called 
the quantum limit, which is the performance achieved in the absence of any 
noise except for the quantum noise inherent to the photodetection prbcess. 
The optimal detection scheme for PPM signals consists of the maximum 
likelihood (ML) detection method. Under the quantum limit, the photo- 
detector output is a Poisson random point process. The probability density 
for the number of detected photons, N t[ts , within an interval [t 1 ,t 2 ] is given 
by [5] 

tj t 2 

Pr[ N tljl =n| = -i-(/X(t)dt)° exp [ - JX( t)dt ]. (1.1) 

t l t. 

The average photon count rate, X(t), is given by 

x « s &ir Pr i *w- 1 1 = tt- w 

where P 0 (t) is the power of the received optical signal, rj is the quantum 
efficiency with which the received photons are converted into photoelec- 
trons, and hf is the photon energy. The sample function density for photon 
absorption times of the Poisson photon counting process over a PPM word 
time is [5] 

or, 

q ~ f X(t)dt 

p[wijw 2 » • • • ^o,qt s =q i = [n x M e ° 

1=1 
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where u x represents the time at which the i th detected photon is registered. 
The likelihood function is defined as the logarithm of (1.3) and can be writ- 
ten as 

QT, N OTj =n 

l=- f \(a)da + £ lnX(w i ) . (1.4) 

0 i=l 


The average photon count rate for PPM signaling is piecewise constant, i.e. 


X(t) = 


{ 


X 3 + x 0 


in the slot containing the PPM pulse 
in all other slots 


(1.5) 


where X s , X 0 are the photon count rates due to signal and background radia- 
tion, respectively. The maximization of (1.4) reduces to 

max {ln(X 0 + X s )N(j_ 1)Tf jX> + § ln(X 0 )N( i _ 1 ) XjjiX| ? 1<J<Q }• (1-6) 


In other words, ML detection of the PPM signal is equivalent to finding the 
time slot which contains the largest number of detected photons. If the 
number of detected signal photons per bit is denoted as n bit , the number of 
detected signal photons per PPM pulse is n bjt Xlog 2 Q. The PPM word error 
probability (WEP) under the quantum limit can be written as 

WEP = exp [ - n bit Xlog 2 Q ]. (1.7) 


The receiver bit error rate (BER) is related to the PPM word error probabil- 
ity by [12] 

BER = 2f^rr WEP - (1,8) 

Equation (1.7) and (1.8) also show that the number of detected signal pho- 
tons per bit required to achieve a receiver BER decreases monotonically as 
the alphabet size of the PPM signal, Q, increases. 
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In a practical optical receiver, the thermal noise and the limited 
bandwidth of the preamplifier have to be considered. The actual output 
from the photodetector can be modeled as the sum of a filtered Poisson ran- 
dom point process and a continuous thermal noise process. The exact max- 
imum likelihood function becomes very complicated. The ML detection 
procedure derived for the quantum limited case may be applied as an 
approximation. The approximate ML PPM signal detection procedure 
which can be implemented consists of integrating the photodetector output 
current over each time slot and then determining which time slot contained 
the largest integrated electron charge. 

There are two types of photodetectors which are suitable for 
wavelengths emitted by AlGaAs lasers (X=780~840nm), silicon photo- 
diodes and silicon avalanche photodiodes (APD). Both of them have quan- 
tum efficiencies of about 80% near a wavelength of 800nm. The former is 
stable and simple to use, but, when the received optical signal power is low, 
the output photocurrent may be buried completely in the thermal noise of 
the photodetector load resistor at the input of the preamplifier. An APD 
photodetector has internal multiplication gain for the primary photocurrent 
as in a photomultiplier tube, and its output photocurrent is much larger 
than that of a ordinary silicon photodiode. Personick [13] has shown that 
use of an APD photodetector can improve the performance of a direct 
detection optical communication system by 10-15 dB in minimum optical 
power required to achieve a BER of 10 -2 -10 -12 . However, the gain of an 
APD is random, which introduces so called excess noise in the output pho- 
tocurrent. A stable and high voltage bias supply is required and some kind 
of temperature compensation may be necessary in order to keep the average 
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APD gain constant. 

The statistics of the APD output photocurrent is rather complicated as 
shown in [14, 15]. In many applications, the APD output is approximated 
as a continuous Gaussian random variable when computing receiver perfor- 
mance. Although the true distribution differs significantly from the Gaus- 
sian distribution at the tails of the distribution functions, it is shown [16] 
that the use of the Gaussian approximation results in only insignificant 
errors when both the background radiation and received signal levels are 
high. It was not known, however, if the Gaussian approximation could still 
be used when the background radiation level is very low. Low background 
radiation may be encountered in a free space optical communication link 
between two satellites when neither the Sun, the Earth, or the Moon is in 
the field of view of the receiver optics [17]. Furthermore, the average 
number of detected background noise photons per PPM slot time decreases 
as the data rate increases, because the intensity of the background radiation 
is constant. Therefore, an accurate and efficient model was needed to 
characterize the APD output photocurrent under conditions of a very small 
number of detected background noise photons per PPM time slot. 

One of the basic problems inherent in using PPM signal formats is 
receiver timing recovery and synchronization. The receiver has to reestab- 
lish not only PPM slot synchronization but also the PPM word synchroniza- 
tion. Offset or jitters in PPM time slot boundaries will cause the ML detec- 
tor to integrate over a time interval which overlaps two PPM time slots. As 
a result, part of the received signal energies is spread into the adjacent 
time slot and hence acts as background noise. Errors in PPM word boun- 
daries will cause the ML detector to compare the integrated APD outputs 
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from time slots that do not belong to the same PPM word. This will result 
in catastrophic detection errors in the regenerated PPM sequence and the 
demodulated binary data sequence. Timing recovery for PPM signaling can 
be established using similar methods for bit and word/frame timing 
recoveries found in conventional digital communication systems that use 
OOK signal formats. Some studies of timing recovery in optical PPM com- 
munication systems have already been done previously 

[18, 19, 20, 21, 22, 23]. However, it is still of interest to search for new 
schemes which are more suitable for PPM signal formats, and simpler to 
implement. The performance of the timing recovery schemes should also be 
well characterized. 

An efficient algorithm was developed that used the nearly exact Webb’s 
approximation [24] for the number of photoelectrons output by an APD. A 
set of upper bounds were found for the resultant computation errors of the 
numerical procedure. It was shown that the results of the Gaussian approxi- 
mation diverged significantly from the results obtained using Webb’s 
approximation when the average numbers of detected background noise 
photons per PPM slot were less than one. 

A laboratory optical communication system was constructed and its 
performance carefully measured. The system used Q=4 PPM signaling at a 
source date rate of 25 megabits per second. A_n AlGaAs laser diode (Hitachi 
HL8314E, X=833nm) was used as the optical transmitter and a silicon APD 
(RCA 30902S) was used as the photodetector. The receiver BER was meas- 
ured as a function of the average number of detected photons per informa- 
tion bit. It was found that use of the Gaussian approximation for the .APD 
output photocurrent was inappropriate when the background radiation level 
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was less than one detected noise photon per PPM slot time. For example, 
when the background radiation level was 0.0122/slot, use of the Gaussian 
approximation gave an optimal average APD gain that was smaller than the 
measured value by a factor of two and a BER that was larger than the 
measured value by 1 to 3 orders of magnitude. The measured values of the 
optimal average APD gain and the receiver BER agreed well with those 
predicted when using Webb’s approximation. 

PPM slot timing recovery in optical communication systems' was 
reviewed and a new scheme was introduced. This scheme, though popular 
in conventional digital communications, had not been studied for optical 
PPM communications. It consisted of a transition detector and a phase 
lock loop (PLL). The transition detector was made of a high speed com- 
parator and a pulse shaper. The output of the pulse shaper consisted of a 
series of short pulses each of which corresponded to the detection of a posi- 
tive transition in the received PPM light pulse train. The power spectrum 
of the transition detector output contained a strong discrete component at 
the PPM slot clock frequency which was tracked by a PLL. The output of 
the voltage controlled oscillator of the PLL then formed a continuous and 
synchronous PPM slot clock. A complete analysis of the phase error in the 
recovered PPM slot clock was presented and it agreed well with the experi- 
mentally measured values. It was shown that use of the recovered slot 
clock at the receiver caused no measurable degradation in receiver sensi- 
tivity if the rms phase error in the recovered PPM slot clock was 1% of a 
slot clock cycle, or the loop noise bandwidth of the PLL was below 10 -4 of 
the slot clock frequency. 
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A novel PPM word timing recovery scheme was invented that used a 
PLL to track the occurrence of back-to-back PPM pulses (PPM word ”Q” 
followed by an ”1”) in the received random data sequences. This PPM 
word synchronization scheme was efficient because there was no prescribed 
synchronization pattern inserted into random data and hence no channel 
use overhead associated with insertion of synchronization patterns. The 
measured rms phase error of the recovered word clock was close to the 
values predicted by the theoretical analysis. A special PLL was devised 
which phase locked the recovered PPM word clock to the slot clock, and 
consequently, corrected phase errors of the word clock in reference to the 
slot clock. Use of the recovered slot clock and word clock caused no 
measurable penalty in receiver sensitivity when compared with a receiver 
that used common clocks at the transmitter and the receiver. This syn- 
chronization system could acquire and maintain the slot and word syn- 
chronization for an received optical signal level as low as 15 average 
detected photons per bit, or 30 photons per PPM pulse. 

The completely self synchronized receiver achieved a BER of 10 -6 at an 
average received optical signal power of less than -63.5 dBm, which 
corresponds to 55 detected photons per information bit. This is believed to 
be the most sensitive receiver yet reported for this type of direct detection 
optical communication system [25]. The receiver sensitivity achieved was 
only 9.1 dB from the quantum limited performance (6.7 photons per bit at 
BER=10 -6 ). The best performance reported so far for coherent detection 
systems that used semiconductor laser was 34 photons per bit using a far 
more complicated phase shift keying (PSK) homodyne detection scheme 
[26], 
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The system described in this dissertation has the potential to provide a 
communication link between two geostationary satellites. The laser diode 
used had an average power limit of 30 mW (14.8 dBm). The total allowable 
losses between the transmitter and the receiver was therefore 78.3 cLB. 
When 20 cm diameters transmitter and receiver telescopes are used, the 
diffraction limited propagation losses between two geostationary satellites 
120° or 72,000 km apart is 66.5 dB. Therefore, a total of 11.8 dB is left to 
account for the optic component losses and the system design margin. An 
even larger margin may be obtained by incoherently combining the laser 
beams from a number of laser diodes to increase the transmitter power, as 
shown in [27, 28]. 

The rest of this dissertation is organized as follows. Chapter II gives a 
detailed analysis of the system performance using the widely used Gaussian 
approximation and the nearly exact Webb’s approximation for the APD 
output photoelectrons. Chapter IE describes the experimental setup includ- 
ing optics configuration, major components, and electronic circuits. The 
procedure and the results of the measurements are also presented. PPM 
slot timing recovery is covered in Chapter IV. Chapter V describes PPM 
word timing recovery. Chapter VI consists of the conclusions and the final 
remarks. 
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Chapter 2. Performance of Optical PPM Receiver with APD Pho- 
todetectors 

In this chapter, the structure and the statistics of the output photo- 
current of an APD photodetector are described. An efficient algorithm is 
developed to numerically evaluate the receiver performance using the nearly 
exact Webb’s approximation for the true distribution function of the APD 
output. Comparisons are made between the use of Webb’s approximation 
and the commonly used Gaussian approximation. It is shown that use of 
the Gaussian approximation seriously underestimates the optimal average 
APD gain and overestimates the bit error probability when the background 
radiation level and the APD leakage current are very low. Perfect timing 
recovery at the receiver is assumed in the analysis. The issue of timing 
recovery and the effects of jitter in the recovered timing signal will be 
addressed in detail in Chapter 4. Chapter 5. 

2.1. APD Photodetectors 

The structure of a silicon APD is shown in Figure 2.1 [29]. A thin 
layer of heavily doped n type (n + ) material is on the top. Next to it is a 
thin layer of less heavily doped p material and a thick layer of very lightly 
doped p material. This second p layer is labeled i for intrinsic because the 
material is so lightly doped that it can be considered to be intrinsic. A 
heavily doped p layer (p + ) is added at the bottom to reduce the contact 
resistance between the device and the outside circuit. 

When a large reverse bias voltage is applied as shown in the diagram, 
depletion occurs at the p-n junction. Because the p layer adjacent to the 
n + layer is thin and less heavily doped, the depletion region can extend 
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through the entire intrinsic layer if the bias voltage is high enough. Due to 
the differences in the doping levels, a large electric field appears across the 
thin junction between the n + and the p layers and a much smaller electric 
field appears across the thick intrinsic layer. Figure 2.1 also shows the dis- 
tribution of the electrical field versus distance inside the APD. 

Incident photons are mainly absorbed in the depleted intrinsic layer. 
Each absorbed photon generates an electron-hole pair. Under the influence 
of the electric field, electrons travel into the high field region where they 
can acquire sufficient energies to excite new electron-hole pairs through 
impact ionization. These secondary electrons and holes may collide with 
other molecules and excite even more secondary electron-hole pairs. The 
process goes on until the electrons and holes leave the high field region. As a 
result, the input primary photoelectron current is multiplied at the output 
of the APD. The multiplication gain is a random number because the 
impact ionizations occur at random. The thickness and the cross section 
area of the intrinsic layer is determined by a tradeoff between the quantum 
efficiency and the response speed of the device. Silicon APD’s can have 
internal quantum efficiencies of nearly 100% at the wavelengths emitted by 
AlGaAs laser diodes (780-840nm) [29]. 

The number of photons which are absorbed in the APD is random and 
the probability that n photons are absorbed from an incident optical field of 
known intensity, P 0 (t) watts, over an interval of [t,t-f-T] is given by a Pois- 
son distribution. The mean number of absorbed photons is given as 
t+T 

n=(r;/hf) J P 0 (t' )dt' , where r\ is the quantum efficiency, and hf is the 

t 

incident photon energy. The probability that n absorbed primary photons 
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give rise to m photoelectrons at the output of the APD was given by McIn- 
tyre [14] and Conradi [15] as, 


nT(- 


P( m | n) 


m 

nr 


eff 


■1) 


k pff m 

m(m-n)! r(— ■ — 4-1+n) 


1— k, 


eff 


[ l+WG-l) | ^ [ (l-k„XG-l) r _ n 


G 


G 


( 2 . 1 . 1 ) 


where G is the average APD gain, k eS =a//? with a and /? the ionization 
coefficients of holes and electrons (i.e. a=probability of ionization per unit 
distance caused by a hole). The probability that m photoelectrons are out- 
put by the APD in response to n average absorbed primary photons can be 
written as 

P(m— 0 | n) = e~ a , (2.1.2) 

P(m | n) = Y, P( m ! n )-^7~ e ~ n > m>l. (2.1.3) 

n— 1 

The noise due to the randomness of the APD gain is called excess noise and 
is characterized by an excess noise factor defined as F=E{m 2 }/E 2 {m}. It 
can be shown that [15] 

F = k e gG+(2-l/G)(l-k eff ) . (2.1.4) 

The excess noise factor increases as the ionization coefficient ratio and/or 

the average APD gain increase and F— +2 as G— »-oo and k eff — >0. Silicon 

APD’s can have values of k e g as low as 0.006 [30]. The average and mean 

square value of the number of photoelectrons output by an APD are given 

by 

nP t+T 

E{m} = G5 = |jf- / P„(t' )dt< 


(2.1.5) 



E{m 2 } = G 2 Fh = 


hf 


/ P 0 (t' )dt' . 

t 


( 2 . 1 . 6 ) 


2.2. Effect of Thermal Noise 

The photocurrent output by the APD need to be amplified for further 
processing. A transimpedance amplifier is usually used to convert the high 
impedance current signal to a low impedance voltage signal. In a practical 
system, there always exist thermal noises from the load resistance seen by 
the APD (input impedance of the amplifier) and the thermal noise of the 
amplifier. These noises can be approximated to appear additive and have a 
Gaussian distribution with zero mean and a constant (white) power spectral 
density within a certain bandwidth. One can model the amplifier thermal 
noise as an equivalent current noise source at the input of a noiseless 
amplifier and use the total input current of the amplifier to analysis the 
receiver performance. The one sided spectral density of the load resistor 
thermal noise current is given by [31] 


d<i R 2 > 4KT F 


, A 2 /Hz 


( 2 . 2 . 1 ) 


df R 

where the bracket means average value, K is Boltzmann’s constant, and T R 
is the absolute temperature. The amplifier thermal noise may be charac- 
terized approximately by an equivalent noise temperature defined as 
T e /T R =<i t £ ta |>/<i R >, where i^j represents the total noise current from 
the APD load resistor and the subsequent amplifier. The spectral density of 
this total noise current may still be described by (2.2.1) with the absolute 
temperature T R being replaced by T e . When the data rates are high, for 
example, greater than 100 Mbits/sec, the bandwidth of the amplifiers has to 
be large and the noise spectrum within the bandwidth can not be assumed 
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to be white. Smith and Personick have given a more accurate model under 
these circumstances [32]. 

As mentioned in Chapter 1, the detection of PPM signals consists of a 
device which integrates the APD output current over each time slot and 
then determines the largest one through comparisons. The output of the 
integrator contains both APD photoelectrons and thermal noise electrons 
accumulated over the time slot. The electron charge due to the total ther- 
mal noise, x t , is a Gaussian random variable with zero mean and variance 
given by 

var{xj = — BT 3 2 = T s , Coulomb 2 (2.2.2) 

where B=l/2T 3 is the noise bandwidth of an ideal integrator, and T s is the 
PPM slot time. 

The advantage of using an APD is that it provides internal gain for the 
primary signal photocurrent which would otherwise be obscured by the 
thermal noise of the load resistor and amplifier. Since the excess noise also 
increases with the APD gain, according to (2.1.4), there exists an optimal 
value for the average APD gain under which the system performance is 
optimal. In the absence of any thermal noise, the optimal APD gain is 
unity and the receiver reaches its ultimate performance which is the quan- 
tum limit. 

2.3. Gaussian Approximation for the APD Output Photoelectrons 

Since the total APD output consists of the discrete number of pho- 
toelectrons described by (2.1.1) and the additive continuous Gaussian ther- 
mal noise, the explicit expression of the probability density function is gen- 
erally too complicated to use in practice. A widely used approximation is to 
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assume the discrete number of photoelectrons output by the APD as a con- 
tinuous Gaussian random variable with the same mean and variance as 
given by (2.1.5) and (2.1.6). Since the amplifier thermal noise is Gaussian 
and independent of the APD output photoelectrons, the total electron 
charge integrated over a PPM time slot, x, is also a Gaussian random vari- 
able with mean and variance given as 


p(x| ^ = ^ e (2 - 3a) 

x = eGn Coulomb (2.3.2) 

2KT T 

a 2 = e 2 G 2 Fn + eI a T s + — Coulomb 2 (2.3.3) 

R 

where n is the average number of absorbed photons in a time slot of dura- 
tion T 3 , and I 3 is the APD surface leakage current. The APD surface leak- 
age current can be treated as a constant dc current because it is the part of 
the APD dark current which dose not get multiplied by the APD gain. 
There is also bulk leakage current in an APD which does get multiplied and 
it can be treated equivalently as a source of background radiation. If one 
defines the signal to noise ratio (SNR) as 


i 2 

SNR = At = 


(eGn f 


a 2 „ „ 2KT e T s 

e 2 G 2 Fn + eI.T, + 


(2.3.4) 


3 s ' R 

The SNR increases as G increases from unity because the numerator of 
(2.3.4) increases while the main contribution to the denominator is still the 
thermal noise, 2KT e T s /R, which is independent of the APD gain. However, 
the excess noise term, e 2 G 2 Fn, increases faster than the the numerator of 
(2.3.4) because F increases with G as shown in (2.1.4). As a result, the SNR 
starts to decreases as G becomes greater than a certain value when the 
excess noise becomes much larger than the thermal noise. In other words, 
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there exists an optimal average AF’D gain that maximum the SNR. The 
receiver performance does not solely depend on the SNR of the input signal 
and hence the optimal value of G that optimize the receiver performance 
may be slightly different. 

The PPM word error probability (WEP) can be written as 

oo x 

WEP = 1 - / p(x | n t ) [ / p(x' | n 0 )dx / ] Q-1 dx (2.3.5) 

-OO -00 

where nj and n 0 represent the average number of detected (absorbed) pho- 
tons in the time slot where a transmitted PPM pulse is present and absent, 
respectively. In a real optical communication system, 

no= a e n s + n bl + n bg (2.3.6) 

and 

ni=n 3 + n b , + n bg (2.3.7) 

where a e is the extinction ratio of the laser diode emitting power between 

ON and OFF states, n s is the average number of detected signal photons in 
a PPM pulse, n bi is the equivalent average number of detected noise pho- 
tons resulting from the APD bulk leakage current within T s seconds, and 
n bg is the average number of detected noise photons within T s due to back- 
ground radiation. The corresponding bit error rate is given by [12] 

BER = 2(^Tj' WEP • (2 ' 3 ' 8) 

The above Gaussian approximation method has been used by Abshire [33 
and Sorensen 34] to compute the receiver performance of optical PPM com- 
munication systems that used APD receivers. 

2.4. Webb’s Approximation for the APD Output Photoelectrons 

A much more accurate approximation for the probability of the 
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number of photoelectrons output by an APD, (2.1.1), was given by Webb 
[24], as 


p (m | n) = 


(2toG 2 F) 1 /2 m-Gn ]3/2 

1 nGF/(F-l) 1 

(m-Gn) 2 


•exp [ 


2nG 2 F(l+ - - m Gn ) 


- ] (2.4.1) 


GnF/(F-l) 

for m values greater than n. Equation (2.4.1) approaches a Gaussian den- 
sity with mean Gn and variance FG 2 n when | (m-Gn) | «Gn. Condi- 
tioned on the number of secondary photoelectrons, the probability density 
function of the total APD output can be written as 


P(x | n) = g P( x I m )P( m | n) , (2.4.2) 

m=0 

where P(m | n) is given by (2.4.1), and p(x | m) is the conditional probabil- 
ity density function of the APD output including thermal noise. Because 
the thermal noise is well characterized by the Gaussian distribution, 
p(x | m) can be written as 

(X-Xm ) 2 

p(x I m) = ~/==~ e ** (2.4.3) 

x m = me + I S T S (2.4.4) 

2KT 

0 s = K + - • ( 2 - 4 - 5 ) 

The final BER can be obtained by substituting (2.4.1) through (2.4.5) into 
(2.3.5) and (2.3.8). 


2.5. Comparison between the Gaussian Approximation and Webb’s 
Approximation 

Use of Gaussian approximation has the advantage of analytical simpli- 
city. It has been shown that the calculated results of receiver performance 
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obtained with the Gaussian approximation were quite accurate in situations 
where the bulk leakage current of the APD is of order of nanoamperes and 
the background radiation levels are relatively high, i.e. when they produce 
more than 10 detected noise photons per time slot [16, 35, 36]. However, it 
was found in our experiments that use of the Gaussian approximation 
resulted substantial underestimation of the optimal average APD gain and 
overestimation of the receiver bit error rate when the number of detected 
noise photons was very small, e.g, n 0 <l. 

Figure 2.2 shows the probabilities of number of photoelectrons output 
by an APD using the nearly exact Webb’s approximation, (2.4.1), and a 
Gaussian probability density function with mean Gn and variance G 2 Fn, 
when n=l and n=120, respectively. The two curves are close to each other 
only when m«Gn and they differ substantially at the tails of the distribu- 
tion which form the main contributions to the WEP given by (2.3.4). The 
smaller n is, the poorer the Gaussian approximation becomes. As shown in 
Figure 2.2, when n=l, the Gaussian approximation is far off from that of 
Webb’s approximation and almost half the Gaussian probability density 
function extends to the unphysical region of m<0. 

The poor estimates of the optimal APD gain and receiver BER that 
result from use of the Gaussian approximation may be understood heuristi- 
cally as follows. Under the extreme case when n 0 =0, the APD output from 
time slots where no light pulse was transmitted consists of only the thermal 
noise. The receiver BER may be reduced to the quantum limit by increasing 
the average APD gain because thermal noise is independent of the APD 
gain and the probability that the primary photoeiectrons are obscured by 
the thermal noise approaches zero as the average APD gain increases. That 


P(m | n) 


-24- 



Figure 2.2. Probabilities of the number of pbotoelectrons output by an APD. The solid 
curves correspond to the Webb’s approximation given by (2.4.1) and the dashed curves 
correspond to the Gaussian distributions with the same mean and variance as those given 
by (2.4.1). The APD parameters used were G=500 and k e j=0.010. The average number 
of absorbed photons are n=120 and n=l, as indicated in the graph. 
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implies an optimal APD gain of infinity and an error probability approach- 
ing the quantum limit. On the other hand, when the Gaussian approxima- 
tion is used, the probability that the APD output from the time slots which 
do contain a light pulse falls below a noise discrimination level, x 0 , is given 
by 


*o 


P(x<x 0 ) = 1 


1 


-co (27re 2 G 2 FE 1 ) 1 /2 

Th t* 


-(x-eGnjJ 2 

exp — — — dx 

2e 2 G 2 FE, 


(2.5.1) 


- / 


-00 




dt 


where 


Th 


x 0 -eGn! 


(2.5.2) 


(e 2 G 2 Fn 1 ) 1 / 2 ' 

It can be seen that Th— ►0 as G — >00 for any value of x 0 because F is pro- 
portional to G. Consequently P(x<x 0 ) -*-1/2 as G-*-oo. That is to say that 
the optimal APD gain cannot be infinite and the minimal error probability 
never reduces to the quantum limit. In reality, n 0 >0 but as long as n 0 is 
small, the Gaussian approximation still underestimates the optimal gain and 
overestimates the error probabilities. The smaller n 0 gets, the bigger the 
discrepancies become. 


Very low noise APD’s can have bulk leakage currents of less than 10 -12 
amperes [37], which is equivalent to a background radiation noise photon 
absorption rate of the order 10 6 /second. In some free space optical commun- 
ication systems, background light levels are negligible when neither the Sun, 
the Earth or the Moon is not in the field of view of the receiver optics [17]. 
The extinction ratio of the laser between ON and OFF states can be kept 
small by biasing the laser diode below its threshold level. Therefore, the 
mean number of detected noise photons per PPM time slot, can be well 
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below unity, at useful data transmission rates of 100 megabits/second and 
higher. The Gaussian approximation does not accurately model the APD 
output under these conditions and the more exact Webb’s approximation 
for the APD output, (2.4.1) has to be used. 


2.6. Numerical Computations Using Webb’s Approximation 

In this section, an efficient algorithm is described which uses Webb’s 
approximation for the APD output photoelectrons. 

2.6.1. Computing BER with Webb’s Approximation 

Substituting (2.4.1) through (2.4.5) into (2.3.5), the PPM word error 
probability can be written as 

00 00 

WEP = £ P(m | nj / ^(x,x m ,c7 2 ) (2.6.1) 

m— 0 -oo 

oo X—X / 

•{ l-[ £ p ( m ' I no)Erf( — )] Q_1 } dx 

m' =o ° 

where 


^(x,x m ,<7 2 ) = -jL- exp [ - ■ ] (2.6.2) 

is the Gaussian probability density function, x m and x m / are the mean 
values of the electron charge for given number of secondary electrons m and 
m' , respectively, and 


Erf(u) = = / 


* e 2 dt . 


-CO 


n/2tt 


(2.6.3) 


Since Erf( 


x-x m / 


x-x m / 


) = 1-Erf( ), x m / «x m , and x varies about its 


a a 

mean during the course of the outer integration in (2.6.1), 




Erf( )«1. We can then use a Taylor expansion for the Q-l power 
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and approximate the expression inside the braces of (2.6.1) and rewrite 
WEP as 

WEP = (Q-l) g P(m | 

m=0 -oo 

oo X— X / 

• E P ( m ' I n 0 )Erf( — )dx (2.6.4) 

m' =0 ® 

This is equivalent to the union bound for the word error probability. 

In a numerical evaluation of (2.6.4), the. two infinite sums and the 
infinite integral have to be truncated to finite terms and limits. Equation 

(2.6.4) can be rewritten as 

M _ Bm _ 

WEP = (Q-l) f] P(m | nj / <f>(x,x m a 2 ) 

m=0 A m 

• £ P(m' I n 0 )Erf( --^L)dx + e (2.6.5) 

m' =o a 

where M, M' , A m , and B m are the new limits of the sums and the integral. 
The resultant truncation error, e, consists of three parts corresponding to 
each truncation. A set of bounds for those errors will be given in the next 
subsection, along with a procedure for choosing the appropriate limits of the 
sums and the integral. 

It should be noted that Webb’s approximation is only valid for 
m>n>0 and terms corresponding to m^nj and m' <n 0 have to be treated 
separately. The contributions from the terms for m<n t in the outer sum of 

(2.6.5) are negligible since P(m | nj)— *0 for values of m which are much less 
than the mean, Gn x . On the other hand, the contributions from the terms 
for m' <n 0 in the inner sum of (2.6.5) cannot be neglected when Qo«l. 
The probability that no secondary electrons are output by the APD, 
P(m / =0 | n 0 ), is given by e -n °, the same as the probability that no primary 
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electron is generated based on the Poisson distribution. If Hq approach zero, 
which occurs at high data rates , low background radiation levels, and low 
APD bulk leakage currents, the probability that no secondary electron is 
output by the APD approaches unity. In other words, when n 0 «l, the 
contributions from the terms for m' =0 in (2.6.5) dominate. Consequently, 
if the truncation error is negligible and n 0 <l, equation (2.6.5) should be 
modified as 

M _ Bm _ 

WEP « (Q-l) £ P(m | n x ) / <f>(x,x m a 2 ) 

111=5,+ 1 A. m 

M ' _ X-Xm' , 

• £ P(m' | n 0 )Erf( -=-)dx 

m' =1 ° 

+ (Q-l) £ P(m | n^-e -50 - / </»(x,x m <r 2 )Erf( - ~~^)dx (2.6.6) 

m=5i+l A. m 

where I S T S is the value of x m < at m' =0, and P(m | n x ), P(m' | n 0 ) are 
given by Webb’s approximation, (2.4.1). The contribution from the second 
term in (2.6.6) increases as n 0 decreases. In the extreme case when n 0 =0, 
the second term in (2.6.6) becomes the only contribution to WEP and the 
first term vanishes because 

S P(m I 5i)/Vj»/) E P(m' I 5o)Erf( 

m=n, A. m m ; — 1 

oo 

< £ P(m’ | S„T,) 

m' =1 

= 1 - e - ' 5 ’— >-0, as n 0 — <-0 (2.6.7) 

2.6.2. The Truncation Errors 

From Eq.(2.6.5), the truncation error is 


£ = £! + £o + e 3 


( 2 . 6 . 8 ) 
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oo 


x-x m / 


e i = (Q-i) E p ( m I Q i) 

m=M+l 

oo 

• / <t>(x,x m ,c T 2 ) £ P(m' I n 0 )Erf(— 

-oo m' =0 

M 

e 2 = (Q-l) £ p ( m I n i) 

m=0 


_ _ 

J <j>{x,x m ,a 2 ) £ P(m' | n 0 )Erf(- 

x>Bn] l x^Afg m' 0 

Bm 


(2.6.9) 

•)dx 

( 2 . 6 . 10 ) 

)dx 


x-x m / 
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e 3 = (Q-l) E p ( m I nJ/^Xjn.cr 2 ) 

m=0 A m 


( 2 . 6 . 11 ) 


x-x m / 


• £ P(m' I n 0 )Erf(- 

m' =M' +1 ° 

A set of bounds for e 2 , and e 3 will be given next. 


■)dx . 


x-x m , 


Since probabilities do not exceed unity and Erf( - — )<1 for any 

<7 


value of x, 


M 

e 2 < (Q-l) £ P(m | n x )- / <£(x,x m ,a 2 )dx (2.6.12) 

m— 0 x>B ra , x<A m 

The integrand in (2.6.12) is now a Gaussian density function and 
A m and B m can be chosen such that the integration is a constant for any 
value of m. Then, 


e 2 < (Q-l) / ^(x^o^dx (2.6.13) 

x>B m , 

A m and B m can be set to be symmetric about x m and deviate by an amount 
which is determined by the specified maximum value for e 2 . For example, 
A m =x m -6cr and B m =x m +6cr for e 2 <(Q-l)X10~ 9 . 

Next, since x m and x m / become further apart as m increases, the value 
of the integration in (2.6.9) decreases monotonically as m increases. There- 
fore, 
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e i < (Q-l) E P(m | n t ) 

m=M+l 

00 co X— X / 

• / ^(x,x M+1 ,a 2 ) £ P(m' j n 0 )*Erf( — )dx 

-oo m' =0 ® 

oo _ — 

oo X—X / 

<(Q-1) f ^(x,x u+l ,a 2 ) £ P(m' | 5^)-Erf( )dx 

-co m' =0 a 

Blktfl — 

M X“X / 

~(Q-1) / «^(x,x m+ 1 ,<7 2 ) j] P(m' | n 0 )-Erf( ^-)dx (2.6.14) 

Ami-i m' =0 a 

Hence it is safe to choose M such that the integral given in (2.6.14) is less 

than a specified value. The integral in (2.6.14) is identical to that in (2.6.5) 

when taking m=M+l, so proper values of M can be determined during the 

course of the numerical evaluation of (2.6.5), rather than through trials of 

separate computations of (2.6.14). Every time the index of the outer sum is 

increased, i.e. a new term is computed and accumulated for the outer sum 

of (2.6.5), the values of the integral within this new term is compared with 

a preset limit, for example, 10~ 9 /(Q-1) for e^lO -9 . The process terminates 

when the value of the integral in the next new term becomes less than the 

preset limit. 

Finding the proper value of M' is lengthy. First, since 

Erf( — — )<1 for any value of x, 

o 

M Bra _ „ co 

*3 < (Q-l) E p ( m I n 1 )/^(x,x m ,cr)dx' £ p ( m ' I n o) 

m=0 A m m / =M' +1 

co ___ 

< (Q-l) E P ( m ' I n o)- (2.6.15) 

m' =M' +1 

Substituting (2.4.1) for P(m' | n 0 ) in (2.6.15) and bounding the sum by the 
associated integral gives, 
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(m' -Gn 0 ) 2 


•exp - 


_ 0 , m' -Gn 0 
2 no G 2 F(l+— 0 


r — Idm' (2.6.16) 
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Since, for m' >Gn 0 , 


, m' -Gn 0 , m' -Gn 0 

(1+— — < (1+ -) = — 

GnoF/(F-l) y GSq ; Gn 0 


then, 
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2n 0 G 2 F(l+— )<2GFm' . 


G?oF/(F-l) 


The exponent in (2.6.16) satisfies 
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G5qF/(F-1) 
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F urthermore, 


m'-Gn, |3/2 . m ' -Gn 0 

%,GF/(F-1) n„GF/(F-l)' ’ 


then, 
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£ < Q-1 ! 


(2to 0 G 2 F) 1/2 m> m' -Gn 0 


E 0 GF/(F-1) 


,] 3/2 


m' -2EUG 

' eXp 1 2GF^ |dm ' 
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2GF 

, r M< _fr ° G 

63 < 2 (Q -l )5 oe no/2F ' V { ^ eMu ( 2 - 6 - 18 ) 

The proper value of M' is found by setting the upper limit of the integral 
in (2.6.18) such that e 3 is smaller than a desired value. Following are two 
examples with two sets of the parameters used in the experiments described 
in Chapter 3. 

Example 1: 

no=l, k eff =0.01, G=500, Q=4. For e 3 <10 -9 , 

2GF 

— — — — < 0.067, or M' >105,000 
M' -n 0 G 

Example 2: 

no=0.01, k eff =0.01, G=700, Q=4. For e 3 <10“ 9 , 

2GF 

——4 — < 0.095, orM' > 133,000 
M' -n 0 G 

2.6.3. Numerical Computation Procedure 

The procedure to evaluate the receiver performance, (2.6.6), is as fol- 
lows. First, the limits of the integral, A m and B m , for a fixed value of the 
index of the outer sum, m, are computed. Corresponding to each value of 
the integration variable, x, the inner sum and then the integrand are 
evaluated. The result of the integration is multiplied by P(m j n t ) and the 
product is then accumulated. The process repeats for all values of m from 
zero to M. The computer program used was written in FORTRAN as listed 
in Appendix A. The subroutines for the integration and error function, 
Erf(x), in (2.6.6) were called from the IMSL library. 

Even more computation time can be saved by increasing the summa- 
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tion indices in (2.6.6) by increments larger than unity with subsequent mul- 
tiplication of individual terms by the value of the increments. Special care 
must be taken for small values of m r where P(m' | n 0 ) changes rapidly as 
m' increases. Experience showed that m' should be incremented by unity 
over the range l<m' <100. Trial and error methods indicated that 
differences in the results from trial to trial could be maintained below 1% 
for increments of m and m' (when m' >100) which were of the order of 
several hundred. 

In the numerical computations, the subroutine for the integration was 
set to have a relative error less than 1 %. The limits of the two sums in 
(2.6.6), M and M' , were chosen according to the procedure given in the 
previous subsection. The proper increments for the indices m and m/ 
(when m' >100) of the sums in (2.6.6) were 2,200 and 500 respectively. 
Therefore about 50 values of m were needed to estimate the outer sum and 
about 350 values of m' (including 100 points for 0<m' <99) were needed 
to estimate the inner sum. As a result, It required about two CPU minutes 
to compute one value of WEP on a VAX 8600 computer. 
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Chapter 3. Experimental Setup and Measurements 

A laboratory optical communication system was built with an AlGaAs 
laser diode and a silicon APD. The system used Q=4 PPM signaling and 
operated at a source data rate of 25X10 6 bits per second. The final receiver 
BER was measured as a function of the average number of detected photons 
per information bit and the results were compared with the numerical 
results from the theoretical analysis of Chapter 2. Perfect timing recovery 
at the receiver was assumed by using the same PPM slot clock and the 
word clock as those in the transmitter. The timing recovery system and the 
receiver BER when using the recovered timing signals will be given in 
Chapter 4 and Chapter 5. 

3.1. Optical Setup 

A schematic diagram of the entire optical setup is shown in Figure 3.1. 
The transmitter consisted of an AlGaAs laser diode and a collimating lens. 
Channel losses between the transmitter and the receiver were simulated by 
a set of neutral density attenuators (~60dB). An optical filter was formed 
with a diffraction grating in Littrow configuration [38] that included the 
focal lens next to the attenuators, the collimation lens in front of the grat- 
ing, and the iris at the focal plane. The lens behind the iris was used to 
focus the filtered light beam to the small photosensitive area of the APD. 

The laser diode (Hitachi HL8314) had a wavelength of 833 am and a 
maximum cw power limit of more than 30 mW. The temperature of the 
diode was kept at 20 0 C by the use of a temperature compensated mount. 
The output power vs. injection current and the modulation levels are plot- 
ted in Figure 3.2. The bias level was set at 40 mA, which was below the 
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Figure 3.1. 
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lasing threshold current in order to keep a high extinction ratio of the laser 
diode between ON and OFF states. The peak power of the laser diode was 
set below the cw power limit, as shown in Figure 3.2, to prolong the lifetime 
of the device. 


An optical filter in front of the APD was necessary to cut out back- 
ground radiation and broad band noise emission from the pulsed laser 
diode. A diffraction grating filter was used for the convenience of adjusta- 
bility of both center frequency and bandwidth. A simple interference filter 
with fixed center frequency and bandwidth could be substituted. The angu- 
lar dispersion of the grating is given by [39] 


d9 

dX 


tan#, 


(3.1.1) 


where 9 is the angle of the diffracted beam, X is the wavelength, and 9 { is 
the angle of the incident beam with respect to the normal vector of the 
grating surface, as shown in Figure 3.1. The bandwidth of the light, AX, 
that passes through the iris can be derived from (3.1.1) as 


A * dX A X f 

AX ss —A# = — — 

d 9 tan#; <t> ms 


(3.1.2) 


where f is the distance from the grating to the focal plane, and ^i ris is the 
diameter of the iris. The resolution power of the grating is defined as the 
ratio of the minimum resolvable wavelength difference to the center 
wavelength and is given by 


— - — = (— sin#i) (3.1.3) 

where W is the width of the grating which is illuminated by the incident 
beam. Since AX mim is a measure of how sharply the frequency response of 
the filter cuts off at the limits of the passband, a good optical filter should 
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have AX min «AX. In the system constructed, <£ iris = 2mm, f=1000mm, 
Q\— 15°, and W c=:20mm. Therefore, the bandwidth of the filter was 
AX «7nm, and the minimum resolvable wavelength difference was 
AX m j n ~ 0.067nm. The grating used had 600 lines/mm and was blazed for 
the wavelength of X=l/jm. The transmission efficiency was measured to be 
30 — 40%. With proper blazing of the groves, a diffraction grating filter can 
achieve an efficiency of 80% [39]. 

The APD used was an RCA-C30902S which had extremely low noise 
and low bulk leakage current. This device had a useful diameter of the pho- 
tosensitive surface of 0.5mm and a frequency response up to 1 GHz. Typi- 
cal rise and fall times were 0.5ns. The typical quantum efficiency at X=830 
was 77% [37]. 

No special attention was paid to the transmission efficiency of the opti- 
cal components. In a practical system, however, all the optics should be 
coated for antireflection, the grating filter should be replaced with an 
interference filter to reduce the size of the system, and the elliptical beam 
profile emitted by the laser diode should be corrected to a circular one 
through the use of a prism pair [40]. 

3.2. Transmitter Electronics 

A block diagram of the transmitter circuit is shown in Figure 3.3. The 
binary data source consisted of a binary word generator (Hewlett Packard 
8018A) which was driven at 25 MHz by an external clock provided by the 
PPM modulator. The electronic circuits of both the transmitter and 
receiver were built with Motorola MCI0K series emitter coupled logic (ECL) 
chips on triple layer wire wrapping boards (Augat UVG series). These ECL 
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circuits had rise and fall times of 2~3 ns and a potential operating speed of 
about 250 MHz [41]. The laser diode driver, which had a temperature con- 
trol unit, was provided by the NASA Goddard Space Flight Center. 

The details of the PPM modulator are shown in Figure 3.4. The tim- 
ing base was provided by the crystal oscillator (Greenray ZY7410) and a 
counter. The crystal oscillator generated the PPM slot clock at 50 MHz. 
The output from the least significant bit (LSB) of the counter provided a 
clock at 25MHz for the binary data sequences. An OR gate was used to 
generate a PPM word synchronization signal which was normally ’’high”, 
and became ”low” when the counter state was ”00”. The circuit worked as 
follows. Shift register #1 was used in series shift mode to temporarily store 
the input binary data. The two binary bits in shift register #1 were 
encoded as a PPM word through a l-to-4 address decoder. The resultant 
PPM word was loaded in parallel into shift resistor #2 at the time slot 
when the PPM word synchronization signal was in ’’low” state, and shifted 
out in series afterwards. The one to one correspondence between the binary 
pattern and the PPM word is shown in Figure 3.5. A timing diagram of 
this PPM modulator is shown in Figure 3.6. The bit time of the binary 
data was 40ns and the pulsewidth of the output PPM signal was 20ns. The 
time delay from the input binary data to the output PPM word was equal 
to the time span of two binary bits, 80ns. 

3.3. Receiver Electronics 

A block diagram of the entire receiver electronic circuit is shown in 
Figure 3.7. Descriptions of this circuit will be divided into the following 
subsections, except for the PPM slot and word clock recovery circuits drawn 
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Figure 3.4. PPM modulator 
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Figure 3.5. Binary patterns and the corresponding PPM patterns 
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received binary data 


Figure 3.7. Block diagram of the receiver electronics. 
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in dashed blocks in Figure 3.7. Common transmitter/receiver clocks were 
used in the experiments and measurements described in this chapter. The 
analysis and experiments concerned with the timing recovery will be given 
in Chapters 4 and 5. 

3.3.1. Signal Amplification 

The output photocurrent from the APD was typical on the order of 
one microampere or less and further amplification was necessary for the fol- 
lowing signal processing. The first stage amplifier, or the preamplifier, had 
to be designed very carefully because the receiver performance depended 
mainly on the noise at this stage. 

The APD and the preamplifier enclosed in the dotted line in Figure 3.7 
were in the same package which was electrically shielded. The package was 
made by Analog Modules, Inc.. The preamplifier had a gain of 28.9mV//iA, 
a bandwidth of 45KHz — 440MHz, and an equivalent input noise current of 
4.4pA/ \/Hz according to the data sheet supplied by the manufacture. The 
load resistance seen by the APD was R=1030fi. The APD gain could be 
adjusted by slightly changing the reverse bias voltage. 

A Mini-Circuits ZFL500 power amplifier with 22dB of gain, and a 
Trontech PI50D power amplifier with 27 dB gain, were used to further 
amplify the signal from the preamplifier. The former had a bandwidth of 
50KHz — 500MHz and a noise figure of 5.3dB. The latter had a bandwidth 
of 30kHz — 150MHz and a noise figure of 5dB. The noise figures were the 
ratios of the variance of the equivalent input noise to that of the input 
resistor of the amplifier. The thermal noises from these two amplifier were 
insignificant compared to the noise of the preamplifier. The resistive 
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attenuator following the amplifiers was used for adjustment of the signal 
levels. 

The lowpass filter after the attenuator, Allen Avionics F2971, was used 
to limit the noise acceptance bandwidth of the circuit. The filter was also a 
part of the matched filter to be described in the next subsection. The cutoff 
frequency of the lowpass filter was 97MHz, which was determined primarily 
by the PPM pulsewidth and the design of the matched filter. The two 1- 
to-2 power splitters in Figure 3.7 were Mini-Circuits ZMSC-2-1, which had a 
bandwidth of 100KHz-400MHz. The lower limit of combined bandwidth of 
the entire circuit was lOOKHz, and therefore, the effects of low frequency 
1/f noise could be neglected. 

3.3.2. The Matched Filter 

The integration of the APD output photocurrent over each time slot 
required by the maximum likelyhood PPM receiver design can be realized in 
hardware with the use of a matched filter. When the pulses output by the 
APD are assumed to have a rectangular pulse shape and a pulsewidth of T s 
seconds, the impulse response of a matched filter should have the same rec- 
tangular pulse shape [42]. The output from such a matched filter at time t 
is equal to the integration of the input signal over the interval [t,t+T s |. If 
the filter is sampled every T s second at the PPM slot boundaries, each sam- 
ple is equal to the integration of the APD output photocurrent over the pre- 
vious time slot. 

The matched filter used is shown in Figure 3.8, which had a tapped 
delay line structure as described in [42]. The time delays consisted of 
RG316 coax cables cut to the proper lengths. The power splitter and the 





Figure 3.8. The matched filter 
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power combiner in Figure 3.8 were Mini-Circuit ZMSC-4-1 l-to-4 power 
splitters, which had a bandwidth of 100KHz-200MHz. The lowpass filter, 
shown as the dashed block in Figure 3.8, was actually the same lowpass 
filter shown in Figure 3.7, since the power splitter and the lowpass filter 
were both linear devices and the order in which the signal passes through 
each of them could be interchanged. The filter was necessary to smooth out 
the staircase output from the power combiner of Figure 3.8. Lowpass filters 
of 50MHz, 97MHz, and 200MHz cutoff frequencies were tested to see which 
gave the smallest receiver BER. The one with 97MHz cutoff frequency was 
finally chosen. 

The response of an ideal matched filter to a rectangular input pulse is a 
symmetric triangle waveform of base 2T S . A value at the top of the triangle 
is proportional to the total energy contained in the received optical pulse. 
In practice, the triangular pulse shape will be distorted because there can 
only be finite number of delay taps in a real implementation of the matched 
filter. The triangular pulse shape output from the matched filter of Figure 
3.8 can be derived as follows. A rectangular pulse input to the power 
splitter can be expressed as 

rec(t) = u(t) - u(t-T g ) (3.3.1) 

where u(t) is the unity step function, i.e. u(t)=0 for t<0 and u(t)=l for 

>0. The signal input to the lowpass filter is then 

y(t) = rec(t) -I- rec(t-T s /4) 4- rec(t-2T s /4) + rec(t-3T s /4). (3.3.2) 

Since the circuit is linear, the output of the lowpass filter can be written as 

m(t) = m 0 (t) + m 0 (t-T s /4) + m 0 (t-2T s /4) + m 0 (t-3T s /4) (3.3.3) 

where m 0 (t) is the response of the lowpass filter to a single input rectangular 

pulse. The Fourier transforms of m 0 (t) can be expressed as 
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M 0 (w) = REC(o;)H(a;) (3.3.4) 

where REC( w ), H(u>) are the Fourier transforms of rec(t) and the impulse 

response of the lowpass filter, respectively. The Fourier transform of (3.3.1) 

is given by 


X(u/) = — . 

J oj J <jJ 

Assuming the lowpass filter has a ’’brick wall” frequency response, i.e. 


(3.3.5) 


it follows that 
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(3.3.7) 


The resultant distorted triangle is obtained by substituting (3.3.7) into 
(3.3.3). Figure 3.9 shows a normalized and centered version of this triangle 
waveform, i.e. m 1 (t)=m(t+0.875T s )/4. The dotted line in Figure 3.9 
represents the perfect triangle pulse shape output from an ideal matched 
filter. 


3.3.3. PPM Detection Circuit 

The received PPM pulse sequence was detected and regenerated from 
the noisy signals output by the matched filter with the circuit shown in 
figure 3.10. The signal output from the matched filter was first split into 
four branches by a power splitter (Mini-Circuit ZMSC-4-1) and then 
delayed by 0, 1, 2, or 3 time slots. Those four signals were sampled simul- 
taneously at the end of the last time slot of a PPM word. Each sample 
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corresponded to the APD output photocurrent integrated over each time 
slot. Six high speed comparators were used to compare every two samples. 
The slot that contained the largest energy was determined from the com- 
parator outputs through the use of the four NOR gates shown in Figure 
3.10. The resultant PPM word pattern was loaded into a shift register at 
the end of each PPM word and shifted out serially afterwards. The sam- 
pling of the matched filter output was done by controlling the time at 
which the shift register of Figure 3.10 was loaded. The sampling could also 
be done by sampling the all six comparator outputs before the shift register 
was loaded. Experiments showed the two method gave the same receiver 
performance. The comparators in Figure 3.10 were Motorola MC1650 one 
bit A/D converters. The time delay between the detection of the received 
noisy PPM signal and the appearance of the regenerated digital PPM 
sequence was one PPM word time, 80 ns. 

3.3.4. The PPM Demodulator 

Figure 3.11 shows a schematic diagram of the PPM demodulator. It 
worked as follows. The pulse positions of received PPM words were 
represented by the state of a two bit counter driven by the PPM slot clock. 
The counter was synchronous with the counter of the PPM modulator 
shown in Figure 3.4. The two binary bits corresponding to a received PPM 
word were decoded by loading the two bits of the counter output into a 
shift register upon the arrival of a PPM pulse. The two binary bits were 
then loaded to another shift register at the end of each PPM word and 
shifted out serially. Figure 3.12 shows a timing diagram of this circuit. 
The time delay between the received PPM word and the output binary bits 
was four PPM time slots, or 80ns. Therefore, the total delay from the input 
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source binary data to the received binary data was 12 time slots or 240 ns. 
This did not include delays caused by the optical propagation and the 
delays introduced by the diode driver, the amplifier, and the coax cables 
which were not part of the receiver. 

The synchronization of the counter in Figure 3.11 was established as 
follows. The output of the OR gate in Figure 3.11 was used as the PPM 
word synchronization signal in the demodulator. A XNOR gate was used to 
compare this signal with the PPM word synchronization signal from the 
modulator in the transmitter. If the two signals disagreed, the outputs of 
the XNOR gate and the AND gate became "low”, which blocked the clock 
that drove the counter. The counter state remained unchanged until its 
state coincided with that of the counter in the PPM modulator. 

3.4. Measurements 

Detailed descriptions about the procedures used to measure the receiver 
BER, the received optical signal levels, and other system parameters will be 
presented in this section. The measurement results will be compared with 
the theoretical results of Chapter 2. 

3.4.1 Receiver BER vs. Detected Photons per Bit 

The binary source data consisted of repetitions of a pseudo random 
sequence 2 20 -l bits long. The received binary data sequences were directly 
compared with the properly delayed version of the source binary data 
sequences with the use of a XOR gate, as shown in Figure 3.13. The output 
of the XOR gate was sampled at the binary data clock rate through an 
AND gate to avoid miscounting consecutive bit errors. The AND gate out- 
put a short pulse every time there was a bit error between the received 
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binary data and the source binary data. A universal counter (Fluke 7250A) 
read out the number of bit errors per second or bit error frequency. The 
receiver BER was determined by dividing the bit error frequency by the 
data rate, 25X10 6 bits/second. The counter recorded at least 100 error 
events before an average bit error frequency was computed. 

The signal levels received by the APD are expressed in terms of the 
average number of detected signal photons per transmitted information bit. 
This is given by the ratio of the number of detected photons per PPM 
pulse, n s , to the number of bits contained in a PPM word, log 2 Q, i.e., 

photo ns/bit = — — — (3.4.1) 

hf log 2 Q 

where rj is the quantum efficiency of the APD, hf is the photon energy, and 
P 1? P 0 are the received optical power in Watts when a PPM pulse is present 
and absent, respectively. 

Since the optical power meter used could only measure average power 
and not the peak power of the PPM pulses, Pj was obtained by multiplying 
the average received optical power by the PPM alphabet size, Q=4, when 
the PPM signal was present. The average received optical power was meas- 
ured by removing the APD and substituting the sensor of an optical power 
meter (UDT350 with UDT 260 sensor head). No other changes were made 
to the optics during the measurements. The APD was then repositioned 
and the receiver BER was again determined. The data were abandoned if 
the two average error frequencies had changed by more than 10%. This 
was done to avoid the effects of any drift in received optical power during 
the course of the measurements. The photosensitive area of the sensor head 
was much larger than that of the APD to ensure that all the light incident 
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on the APD was included in the measurement of the received optical power. 
On the other hand, the background radiation power received by the sensor 
head was much larger than that received by the APD which had a much 
smaller active area. Therefore, this background radiation power had to be 
subtracted from the readings of the optical power meter when measuring 
the average signal optical power received by the APD. The accuracy of the 
optical power meter was assured by comparing its readings with those given 
by other independently calibrated optical power meters. 

The background radiation power actually received by the APD, P 0 , was 
so small that it could not be directly measured with the optical power meter 
used to measure P t . Alternatively, P 0 was determined by measuring the 
noise current output from the APD and using the relationship [32] 


. d<1 ^ = 2^G*F 
df hf 


(3.4.2) 


Here <i n 2 > is the variance of the noise current output from the 
preamplifier, and e is the electron charge. The RF noise power output from 
the preamplifier was measured with a RF power meter (HP-435B with 
8282A sensor). The spectral density of the noise power was determined by 
dividing the total noise power by the bandwidth over which the noise power 
was measured. The spectral density of the noise current variance on the 
left hand side of (3.4.2) was determined by dividing the noise power spectral 
density by the square of the gain of the preamplifier. This method assumed 
the power spectral density of the noise was white, which was confirmed by 
direct measurement of the power spectrum of the APD output signal over 
the range dc to 100 MHz. The extinction ratio of the laser diode between 
ON and OFF states was assumed very high and its effect was ignored dur- 
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ing the measurements. 


3.4.2. Measurements of Other System Parameters 


The optimal average APD gain for a fixed received optical power level 
was determined by adjusting the APD reverse bias voltage until the receiver 
BER reached a minimum. The actual value of the APD gain was measured 
as follows. The APD preamplifier output signal was displayed on an oscillo- 
scope, and the average value of the peak height of the output in time slots 
that contained the received light pulse was estimated by inspection. Since 
the gain of the preamplifier was g==28.9mV/^A and the received optical 
power, P x , could be measured, the average APD gain was obtained by 


G = 


Vs 


(3.4.3) 


eT/Pj/hf 

where V p is the average peak amplitude of the APD response in volts. The 
measured APD gain was only approximate because of the way the peak 
amplitudes of the pulses was determined. 

Amplifier thermal noise can be neglected compared with the APD 
excess noise when a strong cw background radiation is present. The excess 


noise factor, F, can be determined by measuring the noise current variance 
and using (3.4.2). The ionization ratio, k e ff, can be derived from (2.1.4) 
when the excess noise factor is known. Measurements made at several 
values of APD gain gave an average value of k e g very close to 0.010. (0.0097 


± 0.0018). 


The equivalent noise temperature of the preamplifier was obtained by 
measuring the noise power spectral density of the preamplifier when the 
APD bias voltage was reduced to nearly zero (G=0) so that only the 
amplifier thermal noise was present. The one sided noise power spectral 
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density, N, is related to the equivalent noise temperature, T e , by 

N = 4KT e /R A 2 /Hz . (3.4.4) 

It was found that T e =1100°K. The corresponding input rms noise current 

was 7.7pA/\/Hz, somewhat larger than value given on the data sheet of the 

preamplifier of 4.4pA/V < IIz. The extra noise was believed to come from the 

following amplifiers, cables, connectors, and etc.. 

The number of detected noise photons due to the APD bulk leakage 
current, I b , is given by n b |=I b T s /e. The APD bulk leakage current was 
obtained from the data supplied by the manufacturer of the APD as fol- 
lows. The spectral density of the noise current of the APD in the dark was 
given as (d<i d 2 ark >/df) _1 / 2 =3X10 _13 A/v // Hz for G=600. The bulk leakage 
current was related to the noise current by d<i d 2 ark >/df— 2e 2 G 2 FI b [32]. It 
was found that I d =0.0978pA or n b |=0.0122. 

The total APD dark current, I d , was 12 nA as given in the data sheet. 
Thus the APD surface leakage current was I s =I d - GI b =11.9nA. The one 
sided noise spectral density due to this surface leakage current at the input 
of the preamplifier was 2eL,=3.8X10 _27 A 2 /Hz. The noise spectral density 
of the preamplifier was 4KT e /R=5.3X 10 _23 A 2 /Hz. Therefore the APD sur- 
face leakage current had little effect on the measurements and system per- 
formance. 

3.4.3. Results and Discussions 

Figure 3.14 shows the results of numerical evaluations of the receiver 
BER, (2.6.6) and (2.3.8), as a function of the average APD gain for the 
received optical signal levels indicated in the graph. The curve with 
n 0 =12.5 corresponds to I b =0.1nA and P 0 =O. The curve with n 0 =0.0122 



Figure 3.14. Receiver BER vs. average APD gain. The solid curve corresponds to the 
resuits of the Webb’s approximation described in section 2.4. The dashed curve 
corresponds to the resuits of the Gaussian approximation described in seccion 2.3. Other 
parameters used were: APD ionization ratio k e j=0.010. APD load resistance R=1030f2. 
equivalent noise temperature of the preamplifier T e =1100‘K, .APD surface leakage 
current I,— 11.9nA, and the PPM slot time T,=20ns. 
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corresponds to our actual experimental conditions when no actual back- 
ground light was present. The results based on the Gaussian approxima- 
tion, (2.3.1) through (2.3.8) are also plotted (dashed curves) for the purpose 
of comparison. The measured optimal APD gain with n 0 =0.0122 was 700, 
which was in good agreement with the numerically computed results as 
shown in Figure 3.14. The Gaussian approximation works well for n s =12.5 
but very poorly for n 0 =O.Ol22. The validity of equation (2.6.6) was further 
tested experimentally under the condition 1^=1, which was set by arbi- 
trarily introducing a background light source. The measured optimal APD 
gain in this case was 580, close to what the numerical evaluation had 
predicted, as shown in Figure 3.14. 

Figures 3.15 through 3.17 are the numerical results of receiver BER as 
a function of the average number of detected photons per information bit, 
under n 0 =12.5, n 0 — 1, and n 0 =0.0122, respectively. Figure 3.15 shows once 
again the appropriateness of Gaussian approximation under n 0 =12.5. The 
small crosses in Figure 3.16 and Figure 3.17 are the experimental results 
under each background noise level, and they are all very close to the numer- 
ical computation results based on equation (2.6.6) and (2.3.8). These results 
substantiate the validity of the numerical evaluation procedure described in 
the previous chapter using the Webb’s approximation. The Gaussian 
approximation gives a very poor estimate for the optimal APD gain and the 
receiver BER under low background radiation level. At ^=0.0122 and 
n s =110, use of the Gaussian model underestimated the optimal APD gain 
by almost a factor of two and overestimated the bit error probability by 
more than two orders of magnitude over what was actually measured in the 
experiments. Figure 3.14 through Figure 3.17 also show that the 
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Figure 3.15. Receiver BER vs. average number of detected photons per information bit. 
The solid curve corresponds to the results of the Webb’s approximation described in sec- 
tion 2.4. The dashed curve corresponds to the results of the Gaussian approximation 
described in section 2.3. The average APD gain used was G=400, the optimal value 
according to Figure 3.14. The background radiation level was no=12.5. Other parame- 
ters used were the same as those in Figure 3.14. 
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Figure 3.16. Same as Figure 3.15. but at a different background radiation level. ff 0 = 1. 
The Gaussian approximation used G=450, and the Webb’s approximation used G— 500. 
which are the optimal values according to Figure 3.14. Other parameters used were the 
same as those in Figure 3.14. 
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Figure 3.17. Same as Figure 3.15, but at a different background radiation level, 
n 0 = 0.0122. The Gaussian approximation used G=450, and the Webb’s approximation 
used G=700, which are the optimal values according to Figure 3.14. Other parameters 
used were the same as those in Figure 3.14. 
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discrepancies between the results of the Gaussian approximation and the 
nearly exact Webb’s approximation get bigger and bigger as the average 
number of background noise counts decreases. It may be concluded that 
Gaussian approximation should not be attempted when n 0 <10. 

This system has achieved a BER of 10 -6 at a received optical signal 
level corresponding to 55 average detected photons per information bit. The 
best performance reported so far for similar systems was> 80 average 
detected photons per information bit under BER<10 -6 [43, 25]. Therefore, 
the system described in this chapter has achieved the highest receiver sensi- 
tivity yet reported for direct detection optical communication systems that 
used semiconductor laser transmitters and APD photodetectors. 
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Chapter 4. PPM Slot Timing Recovery 

It has been assumed in the previous two chapters that perfect syn- 
chronization between the transmitter and receiver was available. The 
establishment of such timing synchronization will be the major topic of the 
rest of this dissertation. PPM slot synchronization and word synchroniza- 
tion will be addressed in this chapter and the next, respectively. 

4.1. Introduction 

One of the basic problems in a communication system that uses PPM 
signaling is to establish timing synchronization between the transmitter and 
the receiver. As mentioned in Chapter 1, the optimal receiver for a direct 
detection PPM optical communication system consists of a device which 
integrates the APD output current over each time slot and compares the Q 
integrated outputs to find the largest. As a result, the receiver requires 
both PPM slot timing recovery and PPM word timing recovery. Any offset 
or jitter in the recovered time slot boundaries directly affects the results of 
the integrations of the APD output over the time slots. Any part of the 
signal energy that spills over into the adjacent time slot acts as background 
noise and the performance of the entire system is degraded. Incorrect PPM 
word boundaries will cause the receiver to compare the integrated APD out- 
puts from time slots that do not belong to the same PPM word. Further- 
more, PPM word timing errors will paralyze the PPM demodulator and 
cause catastrophic errors in the demodulated binary data sequences. 

In practice, slot timing synchronization requires a regenerated clock 
waveform at the receiver at the slot clock frequency which is synchronized 
with the clock at the transmitter. Once the correct slot timing is provided, 
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PPM word synchronization can be established by a properly initialized 
modulo-Q counter driven by the regenerated slot clock. Details about PPM 
word timing recovery will be given in Chapter 5. This chapter addresses 
only the issues of slot clock synchronization for receivers that use APD pho- 
todetectors. 

Generally speaking, all the technologies for bit timing extraction in 
conventional microwave telecommunication systems can be used for slot 
timing recovery in optical PPM communication systems. Modification is 
necessary, however, to cope with the shot noise of photodetectors that 
appears with the additive Gaussian noise of subsequent amplifiers. The use 
of APD’s improves detector sensitivity significantly but also introduces so- 
called excess noise, which is nonadditive and follows the Conradi distribu- 
tion [14, 15]. Several conventional means of timing recovery have been stu- 
died for use in optical communications [44, 45, 46, 18, 19, 47] and some of 
them are widely used in optical fiber telecommunication systems [48]. Free 
space optical communication systems, however, are characterized by 
transmission channels which have very high losses but can be considered as 
dispersion free. For example, the losses under diffraction limited operation 
between two geosynchronous satellites 120° apart (73,000 km) is about 
66dB if 20 cm diameter telescopes are used as the transmitter and receiver 
antennas. Optical fiber communication systems, on the other hand, are 
dispersion limited rather than propagation loss limited. The received signal 
levels are usually much higher than those of free space optical communica- 
tion systems. Furthermore, timing recovery in PPM systems is different 
than in the case of OOK signal format which are widely used in optical 
fiber systems. PPM signaling has a fixed duty cycle or pulse rate regardless 
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of the data transmitted. 

This chapter is organized as follows. The next section gives a brief 
review of previous studies of slot clock timing recovery for direct detection 
optical communication systems. Section 4.3 discusses transition detector 
type clock recovery schemes, which are widely used in conventional digital 
communication systems but have not yet been studied for optical communi- 
cation systems. A theoretical analysis associated with this system is also 
developed. Section 4.4 discusses the effects of slot timing jitter on system 
performance. Section 4.5 describes the experiments and measured perfor- 
mance of this slot timing recovery system. 

4.2. Review of Slot Clock Timing Recovery in Direct Detection 
Optical Communication Systems 

The response of an APD to an incident optical field can be modeled as 
a filtered compound doubly stochastic Poisson process [5] if there is no 
amplifier thermal noise. The basic problem of timing recovery is to estimate 
slot clock waveforms from the times at which photons are detected from the 
received optical pulse train. One optimal procedure is to use maximum 
likelihood (ML) estimation. Under the assumptions that the APD has an 
unlimited frequency response, the APD output can be modeled as a com- 
pound Poisson random point process with intensity function, X(t), which is 
proportional to the intensity of the received optical field. The likelihood 
function to be maximized can then be expressed as [5] 

P[{N*; t 1 <c7<to} | X(t)j 

ts N ‘r*2 N »r ‘ 2 

= exp { -/X(a)da + £ InfXfa)] + £ ln(P(mj)] } (4.2.1) 

i=l i=l 

where p^N^; t x < t 2 } | X(t)j is the sample function density for the 
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occurrence times, the cTj’s, of the photon absorptions, and N ti ^ is the total 
number of photon absorptions within the interval [tj,t 2 ]. The number of 
secondary photoelectrons generated in the avalanche region of the APD in 
response to each photon absorption is denoted as nip the probability, P(mj), 
is given in [49] and does not explicitly contain the a x . 

If the data pattern is known, i.e. the shape of X(t) is known except for 
the time origin, the maximization procedure requires finding the time origin 
and the set of subsequent pulse boundaries within [t 1? t 2 ] which maximize 
(4.2.1). If the slot period is known, the only parameter to be estimated is 
the time origin. When the data pattern is not known, the likelihood func- 
tion is obtained by averaging (4.2.1) over all possible data patterns. Such 
ML estimation schemes usually require too many computation steps to be 
carried out in real time. Georghiades [20] derived an approximate likelihood 
function that reduced the amount of computation significantly and yet gave 
almost as good synchronization performance as that of the exact likelihood 
function under the conditions of high signal to noise ratio. A digital com- 
puter is still necessary, however, to store the photon absorption times, 
evaluate the approximate likelihood functions for different time slot boun- 
dary shifts, and to then find the maximum through comparisons. In situa- 
tions where the bandwidth of the APD is limited and the data rate is high, 
the amplified APD output current cannot be modeled as a point Poisson 
process, and therefore, the ML synchronization methods described above 
have to be modified accordingly. 

In conventional digital communication systems, a commonly used ML 
synchronization scheme uses an analog device to compute the derivative of 
the likelihood function. The result is used to derive an error signal to 
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correct the frequency of the clock oscillator at the receiver. When the data 
pattern is unknown, an estimated data pattern based on a previously 
estimated timing base may be used to compute the derivative of the likeli- 
hood function. This is referred to as data-added or decision-directed syn- 
chronization [50, 51]. The same method was studied for use in optical com- 
munication systems by Gagliardi [9], It was shown that the computation 
becomes reasonably simple only for a few special pulse shapes, for example, 
pulses of flat top and exponential rise and fall edges. Rectangular shapes, 
which are optimal for peak power limited optical transmitters, such as sem- 
iconductor laser diodes, are not appropriate for this kind of ML estimation 
scheme. 

A popular variant of ML estimation for rectangular pulses is the early- 
late gate scheme [52] though it is not mathematically optimal. The early 
gate integrates the received signal over the first half of the time slot interval 
and the late gate integrates the signal over the second half of the interval, 
based on previously estimated time slot boundaries. If the system is per- 
fectly synchronized, the two values of the integrations will be exactly equal 
on average. Therefore, the difference of the two integrations can be used to 
derive an error signal to adjust the frequency of the slot clock oscillator. A 
detailed study of early-late gate synchronization in an optical communica- 
tion context has been given by Gagliardi [9]. 

A phase lock loop (PLL) can be used to regenerate the slot clock at the 
receiver. Since PPM signals with non return-to-zero (NRZ) rectangular 
pulse shapes do not contain a frequency component at the slot clock funda- 
mental frequency [53], some nonlinearity must be introduced to generate the 
desired frequency component. Several possible nonlinearities have been 
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studied, among them are square law, absolute value, ln(cosh) and fourth 
power [18, 47]. The PLL may be replaced by a much simpler narrowband 
filter. The disadvantage is the resultant static phase error when the input 
slot clock frequency shifts due to Doppler effects or slow frequency drifts of 
the clock oscillator in the transmitter. The static phase error resultant 
from a PLL with an active loop filter is usually negligible. The static phase 
error resultant from a narrowband filter is given by A<^ ~ tan _1 (2Q f Af/f 0 ), 
where Qf is the quality factor of the narrowband filter, Af is the frequency 
shift, and f 0 is the center frequency of the filter. Since the quality factor of 
the narrowband filter is usually high in order to reduce the jitter of the out- 
put slot clock, a small frequency shift of the input slot clock may cause a 
significant static phase error in the output. Andreucci and Mengali [44] 
have studied timing recovery schemes that use narrowband filters for 
return-to-zero (RZ) data. 

Another class of timing recovery schemes, popular in conventional digi- 
tal communication systems, are called transition detectors [52]. They gen- 
erate unipolar output pulses at positive or/ and negative going transitions of 
the input pulses. The power spectrum of the resulting unipolar pulse 
sequence contains a discrete frequency component at the slot clock funda- 
mental frequency, which can then be tracked by a PLL or a narrowband 
filter. This method was used in our system and is demonstrated to be a 
simple, and yet very effective means of establishing receiver slot clock syn- 
chronization. 

4.3. Slot Clock Recovery with a Transition Detector 
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4.3.1. Transition Type Slot Clock Recovery System 

The shot noise of an APD is signal dependent. Unlike the noise in con- 
ventional communication systems, it is neither a stationary nor a Gaussian 

random process. The mean and the variance of the APD output current, 
* 

i A (t), is given by [32] 

°° erjPUr) 

E{i A (t)} - G / + I d ]h A (t-r) dr (4.3.1) 

-oo 

00 67?P ( r) 

var{i A (t)} = eG 2 F / [ — + Idjh^(t-r) dr (4.3.2) 

-00 

where P 0 (t) is the incident optical power in Watts, I d is the total APD dark 
current, and h A (t) is the causal response of the APD to a single photon 
absorption. Other parameters in the above equations are defined the same 
as in Chapter 2. The APD output current has to be further amplified and 
the amplifier noise, which is stationary and Gaussian, has to be added to 
the total noise. As shown in (4.3.1) and (4.3.2), both the signal and the 
noise output from the APD increase with the average APD gain. The 
amplifier thermal noise, on the other hand, is independent of the APD gain. 
Therefore, there exists an optimal value for the average APD gain under 
which the overall receiver performance including the timing recovery is 
optimized. Since the timing error in the recovered clock can always be 
reduced by using a smaller loop bandwidth for the PLL, the average APD 
gain may simply be optimized such that the receiver BER is a minimum, 
assuming perfect slot timing recovery is available. Under such an optimal 
APD gain, the APD shot noise is usually the major contribution to the 
total noise [54]. The amplified APD output signal contains mostly ampli- 
tude noise due to the multiplicative nature of the .APD shot noise. The 
transition detector, by intuition, is the least affected by amplitude 
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fluctuations of all the commonly used nonlinearities that can be employed 
to introduce the desired frequency component at the slot clock fundamental 
frequency in the received PPM waveform. Consequently, it was used rather 
than one of the other methods in our system. 

Figure 4.1 is a block diagram of the transition detector type slot clock 
recovery circuit. The input to this circuit is the amplified APD output 
which consists of rectangular pulses of width T s corrupted by signal depen- 
dent shot noise and lesser amounts of additive Gaussian noise. The lowpass 
filter accounts for the limited bandwidth of the APD and the following 
amplifiers. The comparator acts as a threshold crossing detector which out- 
puts a clean but jittered pulse to the pulse shaper in response to a received 
light pulse. The output of the pulse shaper is a rectangular pulse of width 
T s /2 with its rising edge lined up with that of the input pulse. This choice 
of pulsewidth maximizes the amplitude of the fundamental frequency com- 
ponent at the slot clock frequency [52]. The pulse shaper has a ”dead time” 
so that it cannot be triggered more than once within one time slot. 


4.3.2. Phase Errors of the Recovered Slot Clock 

The theory of PLL tracking systems has been well developed for con- 
ventional communication systems and several studies have been done on 
optical communication systems as well. It has been shown that for a 
sufficiently narrow loop bandwidth, the phase error of the regenerated clock 
waveform has variance given approximately by [18] 

2 ~ 2BlN(u; s ) 

» p (4.3.3)) 

* s 

Here N(u>) is the two sided noise power spectrum, u/ 5 is the angular fre- 
quency of the locked PLL, P s is the power of the sinusoidal component of 
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Figure 4.1. Block diagram of the slot clock recovery circuit. The dashed blocks 
are the same as those in Figure 3.7. 
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the input signal at the lock frequency, and Bl is the loop bandwidth given 
by 

00 

Bi-/ I H PLL (2rf) | 2 df, (4.3.4)) 

0 

where HpLL(2?rf) is the transfer function of the linearized model of the PLL 
given in [52]. 

The phase error of the recovered slot clock obtained with a nar- 
rowband filter can also be approximated by (4.3.3) with the effective loop 
bandwidth Bl given by [52] 

oo 

Bl = / I H nbl [27r(f - f 0 )] | *df, (4.3.5) 

I. 

where H n b f (2^f) is the system function of the filter, and f 0 is the center fre- 
quency. 

Both the signal power and noise spectral density used in (4.3.3) can be 
obtained from the power spectrum of the signal at the input of the PLL. 
Since the received noisy PPM signal is not a stationary random process, the 
power spectrum at the input of the PLL cannot be obtained by taking the 
Fourier transform of its autocorrelation function. The alternative is to use 
the definition given in [53] 

S(w)-Km | X T (w) | ! } (4.3.6) 

T— +oo ^ X 

where X'p(cu) is the Fourier transform of a sample of the signal waveform 
from -T to T given by 

T 

X T (w) = /x(t)e-i wt dt (4.3.7) 

-T 

As in Chen’s approach in [18], the input signal to the PLL is defined as 

x(t) = £ P(t-nQT s — c a T s — r n ) 


(4.3.8) 
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where p(t) represents the rectangular pulse shape, T s is the duration of a 
PPM time slot, c n is an integer from 0 to (Q-l) used to indicate the random 
PPM pulse position, and r n represents random jitter in threshold crossing 
times. If all transmitted data are equally likely, the {c n } form a set of 
independent, identically distributed (i.i.d) random variables. The {r n } can 
also be considered to be i.i.d. random variables. A derivation of the power 
spectrum under these assumptions is given in the next subsection. 

4.3.3. Power Spectrum of the Input to the PLL 

Substitute (4.3.8) into (4.3.7) and let T=NQT S with N an integer, 

X T (w) = S p(o;)e- j ^ nQT * +c “ T ' +r “). (4.3.9) 

n=-N 

Since the c n ’s and the r n ’s are two independent sets of i.i.d. random vari- 
ables, 

E{ | X x (w) | 2 } = | P(w) | 2 £ e jwQT ^ m ' a) E{e jwT ‘ (CnrCB) }E{e J ^ rra ' rB) } 

m,n=-N 

= I P(w) I 2 [ (2N+1) + 

+ | E{e jwT,CB } | 2 | E{e jwr “} | 2 X e jwQT>(m " n) ] 

m^n 

= I PM 1 2 [ (2N+1) (1 - 1 CM 1 2 1 M^M | 2 ) 

+ I CM I 2 I M.M | ! | S I 2 ], (4.3.10) 

n=-N 

where P(a;) is the Fourier transform of the pulse shape p(t), C(u ) and M ^ui) 
are the characteristic functions of the c n and r n , respectively. These are 
given by 

CM = E{e jwT,CB } = 7 T e jwTsI1 (4.3.11) 

^ n=0 


and 
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M»( w ) = E{e Jwr “}. (4.3.12) 

The power spectrum (4.3.6) can be obtained by taking the limit of (4.3.10) 

as T— >co, therefore, 


SH 


p H 1 2 ( i- 1 cm | 2 1 m,m | 2 ) + 


+ I CM I 2 I M.M I 2 lim -L I £ I 2 ], 

T —co 2 1 n== _ N 


N 


(4.3.13) 


00 


00 


Since e_jxn — ^ Yj <5(x-2^k), and <5(ax)=- 


k=-oo 


-5{x), 


1 N 

lim _L I y; e 

2T a t N 


-W* , 2 = Jlm J£- £ ^ra/QTJ] 2 


T-CO 2T 1 QT, n ± 


n=-oo 


47T 2 


00 1 

E lim -^(c^ro/QT,). 


(QT,) 2 T-oo 2T 
Using the relationship that 

6( w )=-L lim 2T— 

27T T— >00 Cc>T 

each term in the sum of (4.3.14) reduces to 


(4.3.14) 


(4.3.15) 


1 .. 11 sin(u>-27m/QT s )T 

Ito ^^M2WQT,)- lim J-. ±2T - -5M2WQT.) 


= -^(^WQTs). 


According to (4.3.11), 


CM I -i| 1-6 


JwT.Q 


J_ i sin ( u,T sQ/ 2 ) 
Q sin(u>T s /2) 


(4.3.16) 


(4.3.17) 


Q ' l - e JwT * 

therefore, | C(27rn/QT S ) | =1 if n=kQ for any integer k, and 
| C(27m/QT S ) | =0 otherwise. Substituting (4.3.14), (4.3.16), and (4.3.17) 
into (4.3.13), the power spectrum of PPM signal becomes 

SM = I p M ! 2 ( q^(i- I CM 1 2 1 | 2 ) + 


+ 


2 * 1 M,M I 2 E «M-2<rk/TJ|. 


(QTJ 


k=-oo 


(4.3.18) 
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The first term in (4.3.18) corresponds to noise and the second term 
corresponds to the discrete frequency components at the slot clock funda- 
mental frequency and its harmonics. Therefore, the noise power spectrum 
can be written as 

N(«) = I P(w) | 2 ( 1 - | C( U ) | 2 | M,H | 2 ). (4.3.19) 

The total power at the slot clock fundamental frequency is given by 

Ps = I MM) I 2 I P(«J I 2 . (4.3.20) 

4.3.4. Characteristic Function of the Threshold Crossing Time 

If the random jitter at the threshold crossing time, the r n ’s, can be con- 
sidered as i.i.d. zero mean Gaussian random variables, then 

1 g 2 ^j2 

MJu) = e 2 ' (4.3.21) 

where a 2 represents the variance of the threshold crossing times which will 

be determined next. 

The input signal to the threshold crossing detector consists of the APD 
output current, i A (t), and amplifier thermal noise, n a (t). The threshold 
crossing time, t c , is the solution to the equation i A (t c ) + n a (t c ) = I th , where 
I ttl represents the threshold level. The APD output current can be 
expressed as the sum of the signal and shot noise, as 

U(t) = s(t) + n s (t) (4.3.22) 

where s(t) = E{i A (t)} and n^t) = i A (t) - E{i A (t)}. At high input optical 

signal levels, the threshold crossing time, t c , varies within a small region 
about its mean value T c =E{t c }. Consequently, it is appropriate to approxi- 
mate s(t) by the first two terms of its Taylor expansion about T c , and write 
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Ith ~ s(T c ) + s' (T c )(t c -T c ) + n s (t c ) + n a (t c ). (4.3.23) 

Since n s (t) and n a (t) are independent and s(T c )=I til , it follows that 

[s' (T C )IV - E{n 3 2 (t c )} + E{n a 2 (t c )}. (4.3.24) 

The variance of the amplifier noise is given as in Chapter 2, 

E{Hs 2 (te)} = — (4.3.25) 
where K is Boltzmann’s constant, T e is the amplifier equivalent noise tem- 
perature, B is the bandwidth of the lowpass filter, and R is the load resis- 
tance seen by the APD. 

The variance of the shot noise, E{n 3 2 (t c )}, is equal to the variance of 
the APD output current, which is a function of t c . Conditioned on t c , the 
APD output current is a filtered compound Poisson process with variance 
given by [5] 

E{a s 2 (t c ) | tj = e 2 G 2 F / X(t)h 2 (t c -t)dt. (4.3.26) 

-oo 

where h(t) is the impulse response of the lowpass filter and X(t) is the pho- 
ton counting intensity function. As far as threshold crossing time is con- 
cerned, one may assume 


X(t) 


_( Xl 

“\*o 


t>t 0 


(4.3.27) 


\Q t<t Q 

where t 0 is the pulse transition time prior to t c . If the laser has a high on- 
off extinction ratio and the background radiation appears to be small, i.e. 


Xi»X 


o> 


tc 


EMtc) I tj « e 2 FG 2 X I /h 2 (t c -t)dt 

*0 

T c -to t c -to 

= e 2 FG 2 Xj [ f h 2 (u)du + / h 2 (u)du]. 
0 T c -to 


(4.3.28) 
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Since t c -T c is very small, the integrand of the second term in (4.3.28) can 
be approximated by 


h 2 (u)«h 2 (T c -t 0 )+2h(T c -t 0 )h' (T c -t 0 )[u - ( T c - t 0 )]. (4.3.29) 

Averaging both sides of (4.3.28) over t c 


“c* 


Te-to 

E{n s 2 (t c )} « e 2 FG 2 X, [ / h 2 (u)du + h(T c -t 0 )h' (T c -t 0 )-a r 2 ]. (4.3.30) 

o 

Next, according to [5], 


s(t) — E{i A (t)} = eG f X(u)h(t-u)du 


-OO 


then 


(4.3.31) 


(4.3.32) 


s' (t) « eGX^^-to) 
where X 1 »X 0 is assumed. Substituting (4.3.25), (4.3.30) and (4.3.32) into 

(4.3.24), the variance a 2 can be solved as 


„ Te " t ° 0 4KT.B 

FG% / h 2 (u)du + 


a T 2 = 


. (4.3.33) 

[GX 1 h(T c -t 0 )] 2 - FG 2 X 1 h(T c -t 0 )h' (T c -t 0 ) 

Equation (4.3.33) may be evaluated by assuming that the lowpass filter 
is ideal and has a ” brick wall” like frequency response, i.e. 

| f | <b 


H(j27Tf) 


-{ 


otherwise’ 


or 


h(t) = 2B* 


sin27rB(t-t d ) 

27rB(t-t d ) 


(4.3.34) 


(4.3.35) 


where t d is the time delay of the filter. 


4.3.5. Average Threshold Crossing Time 

There is an optimal value for the average threshold crossing time T c , at 
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which (4.3.33) is minimized. In practice, it corresponds to an optimal thres- 
hold level applied to the threshold crossing detector. Figure 4.2 is a plot of 
the normalized rms threshold crossing time jitter, a T / T s , as a function of 
normalized threshold crossing time, AT C B, where AT C = T c -(t 0 -f-t d ). In 
other words, AT C is denoted as the time difference between the mean thres- 
hold crossing time, T c , and the midpoint of the rising edge of the input 
pulse, t 0 +t d . The curve in Figure 4.2 is fairly fiat near BAT C = -0.1. In 
fact, <r r /T g increases by only 3% when BAT C changes from its optimal value 
(£=*-0.08) to zero. Consequently, the value BAT c =0 can be used in (4.3.33) 
to a good approximation. This corresponds to a threshold level at one half 
the peak amplitude of the APD output pulses. 

The impulse response of the low pass filter (4.3.35) and its derivative at 
time t=T c -t 0 can be rewritten by substituting the new variable 
AT c =T c -(t 0 +t d ), as, 


h(T c -t 0 ) — 2B- 


sin27rBAT c 

27rBAT„ 


d 27 tBAT c cos27tBAT c - sin2;rBAT c 

-f-h(T c -t 0 = 2B - 

dt 2 ttBAT 2 


(4.3.36) 

(4.3.37) 


F urthermore, 


Te-to Te-t<^ti AT C < 

/ h 2 (u)du = / h 2 (v+t d )dv « / (2B- S - ^ V ) 2 dv 

0 -t 4 -oo 

2xBAT c 

= B + — / (— ) 2 dx, (4.3.38) 

7T o X 

0 . 

where the relationship J ( sin ^ ) 2 dx = ■— has been used. When BAT C — *-0, 

A 


-OO 


x 


equation (4.3.36), (4.3.37), and (4.3.38) reduce to h(T c -t 0 ) = 2B, 

Te-t° 

— h(T c -t 0 ) = 0, and J h 2 (u)du = B. Substituting these into (4.3.33), the 
dt 0 


0.3 




BAT C 

Figure 4.2. Normalized rms threshold crossing ti 
inalized mean threshold crossing time, BAT C . 
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variance of threshold crossing time at BAT c =0 becomes 


a 


2 

r 


FG 2 X! + 


4KT e 

Re 2 


4BGV 


(4.3.39) 


Finally, substituting (4.3.39), (4.3.21), (4.3.20), (4.3.19) and (4.3.11) into 
(4.3.3), the variance of the phase error in the recovered slot clock is given 
by 


<*} = QT s B L -{exp [ 


tt 2 FG 2 a, + 4KT e T s /Re 2 


BT S G 2 (n s ) 2 


- 1 }. (4.3.40) 


where n s «X 1 T s is the average number of detected signal photons per PPM 
pulse. 


The phase error variance, decreases as the bandwidth of the 
lowpass filter, B, increases. However, the effect of false triggering by noise 
spikes is not included in (4.3.40). This may contribute significantly to the 
overall phase error when the bandwidth of the lowpass filter becomes too 
large. In practice, B is also limited by the frequency responses of the APD 
and the amplifiers. In our experiment, B=2/T s was used. 


4.4. Effects of Slot Clock Jitter on the Receiver BER 

Jitter in recovered slot clock at the receiver results in the appearance 
of part of the received signal energy in an adjacent PPM slot, and hence 
acts as a source of background noise. If the integrators are ideal, the 
amount of signal energy which spills over into the adjacent slot is propor- 
tional, on average, to the amount of offset in the slot boundaries. If the 
timing offset of the recovered slot clock is AT S , the fractional offset is given 
by e = AT S /T S , and the effective average number of detected signal pho 
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tons per PPM pulse becomes 

°s' = (l-e)Ss. (4.4.1) 

Assuming the actual number of detected background noise photons is n 0 , 

the effective average number of noise photons in the adjacent slot becomes 

no = n 0 + «n s . (4.4.2) 

The average number of detected noise photons in the remaining Q-2 time 

slots is still n^. The probability of a PPM word error is now given by 
PWE(e) = 1 - 

00 X X 

- fp(x I n s ' +n 0 )[/p(y j n 0 ' )dy][/P(z | n 0 )dz] Q “ 2 dx (4.4.3) 
ooo 

where p(u | XT S ) is the probability density function of the APD output 
given in Section 2.4 for the average number of absorbed photons equal to 
a 3 ' 4-nQ, Oq' , and n 0 , respectively. 

Equation (4.4.3) usually requires excessive computation time when 
evaluated numerically. One alternative is to use the union bound as shown 
below, 

00 00 

PWE(e) « /p(x | n,' +n 0 )/p(y j n 0 ' )dydx 
0 x 

CO CO 

+ (Q-2)/p(x I Hg' +n 0 )/p(y | n 0 )dydx. (4.4.4) 

0 x 

The union bound gives a good approximation when the input signal to noise 
ratio is high and the timing offset is small, i.e. n s /n 0 »l and e«l. A 
detailed description of an efficient numerical evaluation of (4.4.4) is given in 
Chapter 2.. The overall probability of a PPM word error is obtained by 
averaging (4.4.4) over e. If the timing jitter is assumed to have Gaussian 
distribution with zero mean and standard deviation given as <r f = <7^/2 tt, it 
follows that 
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PWE = / PWE(e)- >l e 2<7 ’de. (4.4.5) 

—oo V27T(7 e 

The receiver BER is related to PWE by (2.3.8). 

In reality, none of the high speed integration circuits implemented in 
hardware is ideal and the amount of received signal energy in adjacent slots 
can not be determined simply by using (4.4.1) and (4.4.2). As an example, 
the integrations over the PPM time slots in our experiment were realized by 
sampling the output of the tapped delay line matched filter shown in Figure 
3.8. Ideally, the impulse response of the matched filter should be a rectangle 
of width T s , so that the output in response to a rectangular input pulse is a 
symmetric triangle of base 2T S . A sample at the top of the triangular 
waveform is proportional to the integration of the input signal over the pre- 
vious time interval of length T a . There is no intersymbol interference 
because the triangle waveform vanishes completely at the next sampling 
time. However, due to the limited bandwidth of the system, all three 
corners of the triangle become rounded. If the input rectangular pulse has 
zero rise time and fall time, and the lowpass filter has a ” brick wall” like 
frequency response with cutoff frequency 2/T s , the output of the matched 
filter can be shown as in Figure 3.9. In practice, input rectangular pulses 
have finite rise and fall times, and no lowpass filters have true "brick wall” 
frequency response. In fact, it may be impossible to derive the actual dis- 
torted triangle shape, and consequently determine the exact amount of 
received signal energy that spills over into an adjacent slot due to timing 
offset. However, when the triangle output by the matched filter becomes 
rounded, the values of samples at the three corners of the triangle waveform 
become less sensitive to the sampling time accuracy. Consequently, jitter in 


- 86 - 


the recovered slot clock has less effect on the system performance than 
when a perfect integrator is used. Therefore, (4.4.5) can be used as an 
upper bound for the effect of jitter in the recovered slot clock on the 
receiver performance. The numerical results of (4.4.5) will be presented in 
the next section. To save computation time, only 40 points of the 
integrand were evaluated during the process of integration. 

4.5. Experiments and Performance Measurements 

4.5.1. The Electronic Circuit 

The slot clock recovery system shown in Figure 4.1 was constructed. It 
used a high speed comparator (Motorola MC1650) as the threshold crossing 
detector and a monostable multivibrator (Motorola MC10198) as the pulse 
shaper. The inherent recovery time (dead time) of the multivibrator 
prevented it from triggering more than once within one slot period. Other 
types of pulse shaper may be used also. For example, the circuit shown in 
Figure 4.3 was tested and the resultant phase jitter of the output of the 
PLL was about the same as when using a monostable multivibrator. The 
pulse shaper of Figure 4.3 can work at very high speed although it consists 
of several relatively large size components, such as power splitters and coax 
delay cables. 

The PLL consisted of a double balanced analog mixer (Motorola 
MC 12002), an active second order loop filter, and a voltage controlled cry- 
stal oscillator (VCXO, Greenray N-471D). The active loop filter was built 
with a LM747NC operational amplifier. A detailed circuit diagram of the 
PLL is shown in Figure 4.4 along with a list of the component values. 


c 
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Figure 4.3. Example of another type of pulse shaper. 
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R fb 



R 1 =4.7Kn, R 2 =6.8ICn 
C = 0.10 n F 
R fb = 1MH 

phase detector gain: IQ = 0.159 V/rad. 

VCXO gain: K 0 = -2 ttX 1.63X10 4 rad/sec./V 

natural frequence w n cs* 300 Hz 

damping factor f ^0.6 

loop noise bandwidth Bl — 1000 Hz 


Figure 4.4. The PLL circuit. 
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All the analysis and measurements of the PLL were based on the model 
given by Gardner [52]. The linearized system function of the PLL is given 
as 

Aq(s) _ 2$w n s + iu n 2 
^i( s ) s 2 + 2$J n s + w n 2 

with 0 o (s) the 0j(s) the Laplace transforms of the output and input phase 

processes, f the damping factor, and the natural frequency. These are 
given by w n =(K 0 K d /R 1 C)- 1 / 2 and ? =a; n R 2 C/2, where and K d are the 
gain of the VCXO and the phase detector, respectively. In the system built, 
K^^ttX 1.63 X10 4 radians/second/volt, and K d =0.0159 volt/radian, 
therefore, the nature frequency and the damping factor were 
300 radians /second and f ^0.63. These were verified by the exper- 
imental measurements. The resultant effective loop bandwidth was 
B l ?«1 KHz. The lock-in range of the PLL was measured to be ±40KHz, 
which was, as expected, the same as the tuning range of the VCXO. The 
pull-in range was found to be ±5KHz when a large feedback resistor, shown 
as a dotted line in Figure 4.4, was used to prevent the operational amplifier 
from saturating due to small DC offsets of the input signals. In a real 
satellite-to-satellite communication system, a more sophisticated frequency 
acquisition circuit may have to be added to extend the pull-in range to cope 
with the Doppler effect. 

4.5.2. Phase Error Measurement Apparatus 

Timing errors were measured with the use of an rms phase error detec- 
tor as shown in figure 4.5. A double-balanced frequency mixer (ASM-15) 
was used as the phase error detector. Under the conditions of small phase 
error, the output voltage of the frequency mixer can be approximated as 
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Figure 4.5. Apparatus of rms phase error measurements. 
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Vp d =Kp d sinA0, «K pd A0, where K pd is the phase detector gain, and A0 is 
the small phase difference between the two input signals. The amplified 
output of the phase detector was input to an rms to DC converter (Burr- 
Brown 4341), which had a bandwidth of DC to 450 KHz. The output of the 
phase detector often contained a DC component because it was difficult to 
completely balance out a static phase shift between the two input signals. 
Since the random phase jitter was often very small compared to static 
phase shifts, a coupling capacitor was inserted between the phase detector 
output and the rms to DC converter input. The capacitor and the input 
resistance of the converter formed an RC filter. The capacitor was chosen 
such that the 3dB cutoff frequency was below 0.1 Hz. Therefore, the 
bandwidth of the entire circuit was 0.1Hz — 450KHz. Two identical slot 
clock waveforms with a series of known fixed phase shifts were used to cali- 
brate the circuit. It was found that the circuit had a nearly constant gain 
of 7.62 V/radian for input phase shifts up to ±7t/4 radians. 

Other sources of phase errors that appeared in the phase error meas- 
urements included noise of circuit components, such as the amplifiers, and 
the small phase jitter inherent in the VCXO. These phase errors were not 
included in our model used to compute the phase error, and therefore, they 
should be excluded from the final measurement results. The rms value of 
these phase errors was determined by feeding the jitter free clock signal into 
the PLL and then measuring the phase error of the VCXO output against 
the input clock. Since the phase jitter of the recovered slot clock is statisti- 
cally independent of the VCXO noise and other circuit noise, the following 
relationship holds, 

2 2 t 2 

^ total ’ & 4 > n 


(4.5.2) 
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where cr^tal 1S the total measured rms phase error, <?$ is the rms phase 
error of the recovered clock described by (4.3.20), and cr^ n accounts for the 
phase error due to the noises of the VCXO and the circuit. The jitter of 
the recovered slot clock is then given by 

&<t> — \/ (7 ^ total ~ ~ ^totali 1 ~ ^"( ) 2 ] (4.5.3) 

z ^total 

Measurements revealed that a was about 7.6 X10 -4 radians. Since the 
smallest value of cr^tal measured was about 5X1CT 3 radians, the approxi- 
mation <7^ o ^totai was used throughout the measurements. 

The average APD gain was set to G=700 — 800 which was the optimal 
value for the receiver with a common transmitter/receiver slot clock as 
described in Chapter 2 and 3. The optimal threshold level for the threshold 
crossing detector was set experimentally such that the observed phase error 
achieved a minimum. As predicted in Section 4.3.5, the optimal threshold 
was found to be near one half of the average peak amplitude of the APD 
output pulses, and small deviations from the optimal level resulted in negli- 
gible increases in the measured rms phase error. 

4.5.3. Measurement Results 

Figure 4.6 is a plot of the normalized rms phase error of the recovered 
slot clock, 27T, as a function of the average number of detected signal 
photons per information bit. The small squares represent the measurement 
data and the solid curve was obtained from the theoretical analysis. The 
feedback resistor shown as a dotted line in Figure 4.4 was disconnected dur- 
ing the phase error measurements in order to keep the analysis simple. The 
PLL was brought to lock by momentarily shorting the output and the input 
of the operational amplifier. The rms phase error increased by less than 
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photons/bit 


Figure 4.6. Normalized rms phase error of the recovered slot clock, 
a^/2-n, vs. average number of detected signal photons per bit when 
using a PLL of Bj==lKHz and a narrowband filter of B L =9.5KHz. 
The curves are the theoretical results of Section 4.4. The small 
squares and triangles represent the measurement data. 
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20% when the feedback resistor was connected. The measured phase errors 
of the recovered slot clock using a crystal bandpass filter of effective loop 
bandwidth Bl=9.5 KHz (full bandwidth 19KHz) are also plotted in Figure 
4.6 (small triangles) along with the theoretical results (dashed curve). It is 
shown that the analysis agreed very well with the measurement data. 

Finally, the recovered slot clock was used in the receiver and the sys- 
tem BER was measured as a function of the average number of detected 
signal photons per information bit. This was to show the effect of jitter in 
the recovered slot clock on system performance. The PPM word timing at 
this stage was provided by dividing the recovered slot clock with a 2 bit 
counter. The counter was synchronized as shown in Figure 4.7. The circuit 
worked as follows. The PPM word synchronization signal was a pulse train 
corresponding to counter state ”00”. The two synchronization signals from 
the transmitter counter and the receiver counter were compared with the 
use of an XOR gate. When the normally closed push button switch was 
open, the output of the XOR gate was used to gate the slot clock that 
drove the receiver counter until the two counters were synchronized. When 
the push button switch was released, the output of the XOR gate was over- 
rid and the transmitter PPM word synchronization signal was disconnected 
from the receiver. 

Figure 4.8 shows the results obtained with the slot clock regenerated 
by the PLL (small triangles) along with the results obtained with the com- 
mon transmitter/receiver slot clock (small crosses). The solid curve 
corresponds to the numerical computations using Webb’s approximation as 
described in Chapter 2 which assumed perfect slot timing recovery. The 
dotted curve represents the numerical results for the upper bound given in 


recovered 



Figure 4.7. Initial synchronization of the receiver counter. 
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Figure 4.8. Receiver BER vs. average number of detected signal pho- 
tons per bit when using the slot clock regenerated by a PLL 
(B L =lKHz). The average number of detected background noise pho- 
tons per time slot was n 0 =0.0122. The solid curve is obtained as in 
Chapter 2 assuming perfect slot timing recovery. The dashed curve 
represents the upper bound given by (4.4.5). The small crosses and 
triangles represent the experimental data when using a common 
transmitter/receiver slot clock (perfect slot timing recovery) and the 
recovered slot clock. Other parameters used were T e =llOO°K, 
k eff =0.010, and R=1030fi. 
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(4.4.5) with ( 7 £ determined by (4.3.40). It is shown that the recovered clock 
worked so well that no measurable penalty was observed in terms of 
detected photons per bit at a BER smaller than 10~ 2 . No cycle slipping and 
loss of lock were observed until the received optical power dropped below 15 
photons per bit (30 photons per pulse) which corresponded to a bit error 
probability above 0.01. Figure 4.9 is similar to Figure 4.8 except that the 
bandpass filter (Bl=9.5KHz) was used in place of the PLL. The dotted 
curve is again the upper bound given by (4.4.5) and (4.3.40). Cycle slipping 
in this case occurred only when the detected optical power dropped below 
30 photons per bit (60 photons per pulse). The measured performance was 
shown to be below the upper bound. The limited bandwidth of the 
matched filter made system performance less sensitive to small amounts of 
jitter in the times at which the output waveform of the matched filter were 
sampled. 

Our results for the rms phase errors of the recovered slot clock are 
comparable with the theoretical results obtained by Chen (18] and Ling [19], 
which used different methods but did not consider the effect of the amplifier 
noise and the randomness of the APD gain. For example, under the condi- 
tions of B l T s = 2.0 X10 -5 , 60 average detected photons per information bit, 
and no background radiation, the result in [18] corresponded to 
o^/27r =3.5 X10 -4 using a square law device followed by a PLL. The result 
in [19] with an early-late gate scheme was 0^/277 =8.2 X10 -4 . Our result 
under the same conditions was 2 k =9.8 X10 -4 , and included the effects 
of the random APD gain and amplifier thermal noise. 
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Figure 4.9. Same as Figure 4.8, but with the PLL being replaced by a 
narrowband filter of Bl=9.5KHz. 
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Chapter 5. PPM Word Timing Recovery 

This chapter discusses a novel PPM word timing recovery scheme 
which regenerates the PPM word clock through the detection of back-to- 
back PPM pulse pairs in the received random PPM sequences and the use 
of a PLL. The presence of a recovered FPM slot clock with negligible phase 
jitter is assumed available. 

5.1. Introduction 

PPM word timing recovery is accomplished by establishing a PPM 
word clock which defines the boundaries of PPM words. Generally speaking, 
all the word timing recovery schemes for block codes may be applied to a 
PPM system because PPM signaling is in fact a block code. The most 
popular and easiest scheme to implement requires the inclusion of timing 
markers or synchronization patterns in the transmitted data stream [55]. 
The receiver detects these synchronization patterns and reconstructs the 
PPM word clock. The synchronization patterns usually consist of a 
sequence of code words which is relatively easily distinguished from random 
data. The transmitted data are usually divided into frames each of which 
consists of a synchronization pattern prefix and a number of data words. 
Word synchronization in this case is equivalent to frame synchronization. 
Frame synchronization should be based on a number of detected synchroni- 
zation patterns in order to average out detection errors. 

Several studies have already been done on frame synchronization 
methods for optical PPM communication systems 
[56, 57, 20, 58, 59, 60, 21, 22]. Two suggestions may be made here though 
they are not among the main topics of this chapter. First, synchronization 
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patterns in a PPM system may be constructed such that they are not a 
sequence of legitimate PPM words, e.g. three consecutive pulses followed by 
3(Q-1) empty time slots. Such synchronization patterns are completely dis- 
tinct from the random data sequence, and therefore, can be detected more 
reliably. Second, binary source data, such as computer data, are often 
coded and already contain synchronization patterns in binary form. In this 
case, PPM word synchronization may be combined with frame synchroniza- 
tion of the binary data through the use of the' existing binary synchroniza- 
tion patterns. One can establishs PPM word synchronization by shifting 
the PPM word boundaries slot by slot until the frame synchronization pat- 
terns of the demodulated binary data appear as expected. 

A special property of PPM signaling is that some PPM word subse- 
quences in random PPM pulse trains can serve as synchronization patterns 
to determine the PPM word timing of the entire random PPM sequence 
[23]. One obvious example is a PPM sequence in which a PPM pulse in the 
Q th time slot is followed by a PPM pulse in the first time slot, e.g. the 
second and the third PPM pulses in the PPM sequence of Figure 5.1. This 
back-to-back pulse pair marks the boundary of the two consecutive PPM 
words as well as the boundaries of other PPM words in the sequence 
because there are exactly Q time slots in each PPM word. The probability 
that there exists at least one back-to-back pulse pair in a random PPM 
sequence of finite length is given in [59] although it is quite complicated to 
evaluate. A lower bound for the average frequency of occurrence of back- 
to-back PPM pulse pairs is given by (1/Q) 2 . The probability that a random 
PPM sequence contains at least one pair of back-to-ba^k PPM pulses 
approaches unity as the length of the sequence increases. I I" icfore, PPM 
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Figure 5.1. Example of a self synchronizable PPM sequence. 
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word synchronization can almost always be established in the absence of 
any inserted synchronization patterns as long as the memory of the receiver 
is sufficiently long. Because back-to-back pulses appear at random and 
errors may exist in the regenerated PPM sequence, a memory device and an 
averaging process are necessary in order to regenerate a smooth and accu- 
rate PPM word clock. This synchronization scheme eliminates the channel 
use "over head” associated with the transmission of inserted synchroniza- 
tion patterns. It also eliminates the need for framing the binary source data 
in order to insert the prescribed PPM synchronization patterns. 

When the transmitted data are not equally likely but have some other 
probability distribution, some kind of source encoding may be necessary to 
ensure the appearance of back-to-back PPM pulse pairs. The binary source 
data can be either scrambled so that they appear equally likely, or they can 
be encoded such that the most likely data patterns correspond to PPM 
sequences that contain back-to-back pulse pairs. There are of course other 
PPM patterns besides back-to-back pulse pairs that may be used to recover 
the PPM word clock, for example, two PPM pulses 2(Q-1) time slots apart. 

The rest of this chapter is organized as follows. Next section contains 
a detailed description of this PPM word timing recovery system based on 
the detection of back-to-back PPM pulse pairs and the use of a phase lock 
loop (PLL). A special PLL was also devised that locked the recovered word 
clock with the slot clock and therefore corrected the phase error presented 
in the recovered word clock. Section 5.3 contains the theoretical analysis of 
the random phase errors in the recovered PPM word clock. Section 5.4 
describes the experimental system and the results of actual performance 


measurements. 
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5.2. Word Timing Recovery Based on the Occurrence of Back-to- 
back PPM Pulse Pairs. 

The key components of the system consist of a back-to-back PPM 
pulse filter and a long memory averaging device. Since the PPM sequences 
regenerated from the maximum likelihood detector of Chapter 2 and 3 
assume the existence of PPM word synchronization, they cannot be used 
for PPM word timing recovery. As a result, PPM sequences used for word 
timing recovery have to be regenerated by comparing the received photon 
energy integrated over each time slot against a threshold, as in an OOK 
system. The error probability of these PPM sequences is higher than for 
those of the maximum likelihood detector. Figure 5.2 gives a schematic 
diagram of the circuit. The photodetector output is integrated over each 
time slot with the use of a matched filter. The signal is then compared 
against a threshold level with the use of a comparator. The shift register 
samples the output of the comparator at the end of every time slot, and 
consequently, restores the PPM sequences. The two bit shift register and 
AND gate form a filter which passes only back-to-back pulse pairs. The 
PLL tracks on these pulses and generates a synchronous PPM word clock at 
the frequency f=l/QT s (Hz), with T s the time duration of a PPM time slot. 

The word clock output by the PLL contains random phase errors 
because the occurrences of back-to-back pulse pairs are random. However, 
since the PPM word clock should be in phase with the slot clock, phase 
errors of less than half a slot clock period can always be corrected by 
further phase locking the word clock to the slot clock. Here the recovered 
PPM slot clock is assumed to have much less phase jitter than the 
recovered word clock. Figure 5.3 shows the circuit that performs the phase 
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Figure 5.2. PPM word timing recovery circuit. 
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locking. The circuit is in fact a different kind of PLL with an exclusive 
NOR gate acting as the phase detector, an AND gate as the loop filter, and 
a counter driven by the slot clock as the voltage controlled oscillator. The 
input signal to this circuit consists of the word clock generated by the PLL 
of Figure 5.2. The exclusive NOR gate compares this input word clock with 
a second word clock obtained by dividing the slot clock with the use of the 
-i-Q counter. The two word clocks have the same waveform except for a 
possible time shift which may have resulted ’ from the phase error in the 
input word clock and the randomness of the initial phase of the second 
word clock generated by the counter. The error signal output by the 
exclusive NOR gate is then used to gate the slot clock and subsequently 
shift the initial phase of the resultant word clock until the two are in lock. 
For example, assuming the two word clocks are initially out of phase 7 r radi- 
ans (2T S second time delay), the exclusive NOR gate will output a ’’low” 
state and block the first two ticks of the slot clock that drives the counter. 
The word clock output from the counter is then delayed by two slot clock 
periods and then is locked to the input word clock. The time for this cir- 
cuit to lock up is less than one half a word clock cycle or QT s /2 seconds. 
As shown in the timing diagram Figure 5.4, the error signal output by the 
exclusive NOR gate has no effect on the gated slot clock when the loop is in 
lock and the time shift of the input word clock due to random phase errors 
is within +T s /2 second. In other words, when the loop is in lock, the word 
clock output by the counter is always in phase with the slot clock if the 
phase error of the input word clock is within +tt/Q radians. The input 
word clock may be delayed by T s /4 seconds to change the maximum 
correctable phase error to ±tt/ 2Q rad. 
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Figure 5.3. Phase lock between word clock and slot clock. 
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Figure 5.4. Timing diagram of the circuit in Figure 5.3. 
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5.3. Phase Error Analysis 


5.3.1. Phase Error Variance of the PLL 


The phase error variance of the word clock output by the PLL of Fig- 
ure 5.2 can be computed in the same way as in Chapter 4.3.2. The phase 
error variance, o$, is given by 




2B l N(o/ w ) 


(rad. 2 ), 


(5.3.1) 


where N(cu w ) represents the noise spectral density at the word clock fre- 
quency w w =27r/QT s , and P w represents the power of the sinusoidal com- 
ponent at the fundamental word clock frequency. The loop bandwidth of 
the PLL, B l , is defined by 


00 

B L = I I H PLL (27rf) | 2 df (Hz) (5.3.2) 

o 

where HpLL(27rf) is the transfer function of the linearized model of the PLL 
given in [52]. 

Both N(a; w ) and P w in (5.3.1) can be obtained from the power spectrum 
of the input signal to the PLL of Figure 5.2, defined as 

T 

W(o») = lim -^-E{ | fy(t)e~> ut dt | 2 } (Watts/Hz). (5.3.3) 
T — >co 2 1 

The input signal, y(t), to the PLL consists of rectangular pulses of width T s 
seconds, each of which corresponds to the detection of a back-to-back pulse 
pair. The input signal can be modeled as 


00 

y(t) = E ^(t-nQTs) (volts/ y/R,), (5.3.4) 

n=-oo 

where p(t) is the pulse shape function, 3^=0, 1 is a random variable 
corresponding to the presence or absence of an input back-to-back pulse 
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pair, and Rj is the input resistance of the PLL. When the transmitted ran- 
dom PPM words are mutually independent and equally likely, the probabil- 
ity of the an’s can be written as 

P r ( a n = l ) = (-^-) 2 

Pr(a n =l | a n±1 =l) = 0 . ( 5 - 3 - 5 ) 

A detailed derivation of the power spectrum (5.3.3) is given in the next sub- 
section. 


5.3.2. Power Spectrum of the Input to the PLL 

Taking T=NQT S with N an integer, the Fourier transform of (5.3.4) 
from -T to T becomes 

T 

/yOOe-^dt = P(w) £ a n e- jwnQT * (5.3.6) 

-T n=-N 

where P(u/) is the Fourier transform of p(t). Since E^an}— 1/Q 2 , 
and E{a n a m }=E{a n }E{a m }=l/Q 4 for | n-m | 7 ^ 0 , 1 , 

T 

E{ | / yftje-i-dt | 2 } = | P( W ) | 2 E E{a n O'T i “ QT,( ' wn) (5-3.7) 

-T m,n=-N 

= | P( W ) I 2 [ (2N+l)-i- + E 1 

W, m,n=-N Q 

| n-m | 7^0,1 


= I PH I 2 [ (2N+1)( ^-~)] 

2N / -jwQT, , JwQT,\ i 1 | -jwQT,n i 2 ] 

+ 11 

The power spectrum (5.3.3) becomes 

T 

WM = lim — l=rE{ | /y(t)e-i^dt | 2 } 

T— co 21 T 


(5.3.8) 
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+ ix” 2T I I ■ 2 1 • («••) 

As in Section 4.3.3, the last term of (5.3.9) can be written as 

>>” W 1 £ e ' i " QT '“ i 2 = 7^ 2 ^2A/QT s ). (5.3.10) 

T-*00 z I n== _ N (QTgf k= _oo 

Substituting (5.3.10) into (5.3.9), the power spectrum becomes 


WM = | P(w) | - ^r-^-QTJl 


+ I P(«) 


2tt 


00 


E <5(w-27rk/QT 3 ) 


(5.3.11) 


Q 4 (QT S ) 2 k ±oo 

The first term in (5.3.11) represents the noise spectral density and the 
second term represents discrete spectral components at the PPM word clock 
frequency and its harmonics. Therefore, the power of the sinusoidal signal 
component and the noise spectral density at the frequency u> w — 27r/QT s 
can be expressed as 


p.- |PK)I 2 - 


Q (QTJ 2 


(Watts) 


(5.3.12) 


NK) = I P(« w ) I 2 7 ^p(~y( 1- ~ 2 ")] (Watts/Hz). (5.3.13) 
Substituting (5.3.12) and (5.3.13) into (5.3.1), the phase error variance of 
the regenerated PPM word clock becomes 

<v 2 = QB l T s (Q 2 -3) (rad. 2 ). (5.3.14) 


5.3.3. Design of the PLL 

Since the circuit in Figure 5.3 can correct for phase errors up to ±7r/2Q 
radians, the loop noise bandwidth of the PLL should be designed such that 
00«7r/2Q. The probability of having uncorr^Hble phase errors is there- 
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fore negligibly small. However, Bl should be large enough for the PLL to 
have a reasonable tracking capability. The tracking range of the word tim- 
ing recovery system should be close to that of the PPM slot timing recovery 
system. 

The acquisition time of the entire word clock recovery system is deter- 
mined mainly by the PLL in Figure 5.2. The PLL can lock up within one 
cycle if the frequency of the received word clock and the free running fre- 
quency of the VCXO are within the lock-in range [52]. If the input word 
clock frequency is out of the lock-in range, some kind of added frequency 
acquisition circuit may be necessary and the acquisition time becomes lim- 
ited by the circuit used. Details about added frequency acquisition circuits 
can be found in [52]. The frequency uncertainty of the input PPM word 
clock is caused by the jitter of the clock oscillator in the transmitter and 
frequency shifts due to the Doppler effect. The jitter of a crystal clock 
oscillator is typically 0.005% and the percentage Doppler shift of the PPM 
word clock is equal to the ratio of the relative velocity between the 
transmitter and the receiver to the velocity of light. 

The PLL may slip cycles or even break lock after acquisition is com- 
pleted. The average cycle slipping time may be used as a lower bound for 
the average time during which the receiver can maintain its PPM word syn- 
chronization. It is shown in [52] that the time between cycle slips is 
exponentially distributed with mean given by 

T a v= (1/B l) exp(jrSNR L ), (5.3.15) 

where Bl is the loop bandwidth given by (5.3.2), and SNRl 1s the loop sig- 
nal to noise ratio defined as 
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SNB ^ ' (5 ' 3 ' 16) 

5.3.4. Effects of Detection Errors 

Errors in the PPM sequences output by the comparator cause the sig- 
nal to noise ratio to decrease at the input of the PLL. There are two types 
of errors, namely, misses and false alarms. Misses in the detection cause the 
signal power to decrease and false alarms cause the noise to increase. A 
quantitative analysis is given as follows. 

First, if there are only misses present, (5.3.5) should be rewritten as 

Jpr(%=l) = 

}pr(a,=l I iVi=l) = 0 (5 ' 3 ' 17) 

where p ms is the probability of a miss. A derivation similar to the one in 
Section 5.3.2 shows the signal power and the noise power spectral density 
become 

Pw' = ( 1-Pms) 4 Pw (Watts) (5.3.18) 

Q 2 -3(i-p ) 2 

N' K) = (1- Pms ) 2 J; - NK) (Watts/Hz) (5.3.19) 

Q -3 

where P w and Nfu;^) are the signal power (5.3.12) and noise spectral density 
(5.3.13) in the absence of detection errors. 

When there are both misses and false alarms, the signal power is still 
given by (5.3.18), but the noise spectral density increases because some of 
the detected back-to-back pulse pairs are erroneous. An erroneous pulse 
pair may consist of either two false alarm pulses or a false alarm pulse and 
a correctly detected PPM pulse. The probability that two false alarmed 
pulses appear in a pair of consecutive time slots is given as 
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Pi= 


I 



p,I (-^p-) 2 


if the two slots are in one PPM word 

(5.3.20) 

if the two slots are in two PPM words 


where pf a is the probability of a false alarm. The probability that one false 
alarmed pulse and a correctly detected PPM pulse appear in a pair of con- 
secutive time slots is given as 


P 2 = 


1 

2 q' p “’ 

2 Lh 

Q Q 


Pfa> 


if the two slots are in one PPM word 

(5^3.21) 

if the two slots are in two PPM words. 


The factor of two in (5.3.21) accounts for the two different orders possible 
for the two pulses. If the probability of a false alarm is small, the effects of 
the appearance of incorrect back-to-back pulse pairs from two false alarmed 
pulses can be neglected since Pi«P 2 as Pf a «l. Under these assumptions, 
the extra noise input to the PLL of Figure 5.2 due to false alarms alone can 
be approximated as 


n fa (t) = E bnlP^-CiiTs-aQTs) + b nrP[ t -( c n+ 1 ) T s- n Q T s] (5-3.22) 

n=-oo 

where b nl , are two random variables, equal to zero or one, which indi- 
cate that a false alarmed pulse has appeared to the left or right side of the 
correctly detected PPM pulse, and c n is an integer from 0 to (Q-l) used to 
indicate the pulse position of the n th PPM word. For equally likely ran- 
dom data, the c n ’s form a set of i.i.d. random variables with uniform distri- 
bution. The probabilities of b n | and b nr can be approximated as Pr(b nl =l) 
= Pr(b nr =l) = p fa . 

The power spectral density of this extra noise signal at the PPM word 
clock frequency can be derived as follows. The Fourier transform of (5.3.22) 
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from -T to T is given by 
X 

/ n, a (t)e-i" t dt = P(w) £ e- i "° QT, -e- i “- T -(b nl + b 1 , r e- i " T ') (5.3.23) 

-T n=-N 

where T=NQT S . The power spectrum at frequency a> w =27r/QT 3 can be 
written as 


N fa K)= lim -^rE{ | / n fa (t)e JWwt dt | 2 } 

I — ►OO L X. ^ 


= | Pk) | 2 lim -L £ E^ 2 **'-^} 

T-oo Z i n m= _ N 

•E{(b„,+b„ e -i 2 '/9)(b mi +b mI e’ 2 ’/Q)}. 


For m^n, 


E ^ e -j2 J Kc B -c nl )/Q^ = | E { e -J2irc/Q}. j 2 

= | i WQ | 2 = o . 

^ c B =0 


For m=n, 


E{e -j2 ^ Cn-c,n ^ Q } = 1 


and 


(5.3.24) 


(5.3.25) 


(5.3.26) 


E{(b a ,+b„e-i 2 '/«)(b ml +b mI ei 2lr /' : l)}= 2p, s +2p, 2 cos^ ~ 2p,. . (5.3.27) 

Substituting (5.3.25), (5.3.26) and (5.3.27) into (5.3.24), 

NfcK) - 2 Pll -L- I Pk) I 2 = 2p,»-§^-NK) , (5.3.28) 

QT S Q -3 

where N(o; w ) is given by (5.3.13). The total noise power spectral density at 
the input to the PLL of Figure 5.2 can be considered approximately as the 
sum of (5.3.19) and (5.3.28), i.e. 

N Kl >») = [ (l-p mi ) 2 Q2 -g^ + 2P, 4 ^ ] Nk) , (5.3.20) 
with N(w w ) given by (5.3.13). The phase error variance at the output of the 
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PLL can be obtained by substituting (5.3.18), (5.3.12), (5.3.29), and (5.3.13) 
into (5.3.1), as 



QB l t s [ 


Q 2 

(i-p ms ) 2 


-3 + 


2PfaQ 4 

U-Pms ) 4 


(5.3.30) 


5.3.5. Probabilities of Miss and False Alarm 

As in Chapter 2, the input signal to the comparator can be considered 
as the electrical charge output from the APD and an equivalent amount of 
amplifier thermal noise accumulated over a time slot. The probabilities of a 
miss or false alarm are given by 

x T 

Pms= / P( x I ^i)dx (5.3.31) 

-00 
CO 

Pfa~ / p( x I °oM x (5.3.32) 

XT 

where x^ is the threshold level, n L and n 0 are the average numbers of the 
detected photons over a PPM time slot when a PPM pulse is present and 
absent, respectively, and p(x | n) with 5=?!,^ is the appropriate probabil- 
ity density function. The average symbol error rate of the regenerated 
PPM sequence is related to p ms and p fa by 

Pe= -^-Pms + -^"Pfa ' (5.3.33) 

The average APD gain should be determined by the PPM detection circuit 
rather than the word timing recovery circuit. The threshold levels applied 
to the comparator may be varied such that the symbol error probability, 
(5.3.33), is minimized. 

The probability density function, p(x | n), using Webb’s approximation 
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is given in Chapter 2.4 as 


00 


P(x | n) = £ p(x | m)P(m | 5) 
m=0 

_ (X-Xm) 8 

p(x | m) = — JL e 2ffS 


with x m = me + I S T S and a 2 — (el s + 


(5.3.34) 

(5.3.35) 


2KT e 

R 


■)T S 


and 


P(m=0 | n) = e- ff 
P(m | 5) = 


(5.3.36) 


(2TmG 2 F) 1 / 2 , m-Gn , 3 / 2 
1 nGF/(F-l) J 
(m-GE) 2 


•exp - 


2nG 2 F(l- 


m-Gn 


m>0 


(5.3.37) 


■) 


GnF/(F-l) 

where m represents the number of secondary photoelectrons generated by 
the APD, Ig is the APD surface leakage current, e is the electron charge, K 
is Boltzmann’s constant, T e is the equivalent noise temperature, R is the 
load resistance seen by the APD, F is the excess noise factor given by F= 
k e gG + (2-l/G)(l-k eff ), k eff is the ratio of the APD ionization coefficients of 
electrons and holes, and G is the average APD gain. 

When (5.3.31) and (5.3.32) are evaluated numerically, the infinite sum 
of (5.3.34) can be approximated by the first M terms. The resultant errors 
may by upper bounded as follows. Substituting (5.3.34) through (5.3.37) 
into (5.3.31), 

XT 


00 


£ms= I E P( x I m ) P ( m I n l) dx 

—oo m=M ra 


X T 


00 


< / p(x I M ms )dx- £ P(m I n t ) 

m=M m i 


-oo 
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/ 

-00 


1 ~T 

>/5F 6 


dt . 


(5.3.38) 


Next, substituting (5.3.34) through (5.3.37) into (5.3.32), 


00 _ 
e fa= J E p( x 1 m )P( m | n 0 )dx 

Xt m=M fl 

oo _ 00 

< El P( m I n o) < / p ( m I n o)dm . (5.3.39) 

m=Mf 1 M^-l 

As shown by (2.6.15) through (2.6.18) of Chapter 2, for M fa »n 0 G, ' 


2GF 

Hoe" o/2F N/F Mf ‘‘ 1 ;“ oG , , 

e fa < 7 1 • f — =re _1 / u 


(5.3.40) 


The values of M ms and M fa should be chosen such that e ms «p ms and 
e fa«Pfa based on (5.3.38) and (5.3.40). For example, the values used to 
obtain the numerical results presented in the next section were 
M ms >MT+ 30,000 and M fa > 100,000 for e ms ,e fa <10 -8 . The parameters used 
were: G=800, k eff =0.010, n 0 =0.0122, I s =11.9nA, T r =1000“K, R=1030n, 
and T 3 ==20ns. 

The error bounds derived above can also be used to compute the error 
probability of communication systems that use OOK signaling. 


5.4. Experiments and Measurements 

The circuits shown in Figure 5.2 were constructed. The matched filter 
was the same one shown in Figure 3.7. The comparator used was a 
Motorola MC1650 and the shift register was a Motorola MC10141. The 
PLL was constructed as shown in Figure 5.5, which consisted of a phase 
detector (MC12002), an operational amplifier (LN747NC), and a VCXO. 


-117- 


Rfb 



R^l.OKH, R 2 =10Kn 
C = 0.10 juF 
Rfjj = 4.7Mf2 

phase detector gain: K d == 7.64 X10" 3 V/rad. 

VCXO gain: K 0 = -27TX1.08X10 3 rad/sec. /V 

nature frequence w n ~ 110 Hz 

damping factor f «0.4 

loop noise bandwidth B L cs* 380 Hz 


Figure 5.5. The PLL circuit for PPM word timing recovery. 
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The VCXO (Greenray ZN-123Z) had a center frequency of 50MHz which 
was Q=4 times the word clock frequency. Therefore, a two bit counter 
(MC10136) had to be used to divide the frequency down to the correct value 
of 12.5 MHz. The OR gate in Figure 5.5 served to generate a unsymmetric 
word clock, which was ”low” only when the counter output was ”00”, the 
same as the PPM word synchronization signal defined in Chapter 3. The 
resistors and capacitor values of the loop filter were R^lKH, R 2 =10KO, 
and C=0.1^F. The resultant loop noise bandwidth was B L «380Hz accord- 
ing to [52]. The standard deviation of the phase error of the output word 
clock was (70=0.020 radian according to (5.3.14). 

The circuit of Figure 5.3 was implemented equivalently as shown in 
Figure 5.6. The symmetric input word clock in Figure 5.3 was changed to 
the the PPM word synchronization signal generated by the PLL in Figure 
5.5. The output of the ~Q counter in Figure 5.3 was also changed to the 
PPM word synchronization signal. No separate circuit shown in Figure 5.6 
was actually built since it was already part of the PPM demodulator shown 
in Figure 3.11. Since all phase errors within ±7r/2Q=±0.39 radians could 
be corrected, the final word clock contained virtually no phase errors at all. 

The average symbol error rate of the detected PPM pulse sequences, 
p e , was measured with apparatus similar to that shown in Figure 3.13. The 
results are plotted in Figure 5.7 (small squares) as a function of the average 
number of detected signal photons per information bit. The threshold level 
applied to the comparator was optimized for each received optical signal 
power level such that the measured symbol error rate achieved a minimum. 
The numerical results of (5.3.33) are shown as the solid curve in Figure 5.7. 
At the optimal threshold levels, of the total number of symbol errors in the 
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Figure 5.6. Implemented circuit of Figure 5.3. 
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Figure 5.7. Probability of symbol errors of the PPM sequences at the 
shift register of Figure 5.2 vs. average number of detected photons per 
information bit at optimal threshold levels. The small squares 
represent the measurement data and the solid curve represents the 
theoretical result given by (5.3.33). The parameter values used were: 
background radiation level n^ = 0.0122. average APD gain G=300. 
k eff =0.010. .APD quantum efficiency r?= 77%, .APD load resistance 
R=1030 ohm. equivalent noise temperature T=1100°K. .APD surface 
leakage current I s =11.9nA and PPM slot time T s =20ns. 
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detected PPM pulse sequences, l/Q are misses and (Q-l)/Q are false 
alarms. In other words, we can use Pms^Pe/Q and Pf a ~Pe(Q“l)/Q when 
evaluating (5.3.18), (5.3.29) and (5.3.30). 

Figure 5.8 is a plot of normalized optimal threshold levels as a function 
of the average number of detected photons per bit, where the normalization 
is with respect to the average peak values of the matched filter output. It 
was found that the optimal threshold levels varied from 24% to 32% of the 
peak value of the matched filter output. The average DC level of the 
matched filter output was 25% of the peak value and independent of the 
received signal level. Since the matched filter and the comparator were AC 
coupled, the effective threshold levels were automatically adjusted close to 
the optimal values if the threshold level was fixed at zero volts. 

Figure 5.9 shows a measured power spectrum of the input signal to the 
PLL of Figure 5.2 at a received signal level corresponding to an average of 
55 photons/bit. The false alarm and miss probabilities were negligibly 
small (~10 -4 ) at this received signal level. The measured data agreed well 
with the results obtained from the theoretical analysis (5.3.12) and (5.3.13) 
for rectangular ECL pulses and a 500 input resistance. The calculated 
power of the sinusoidal word clock component was -22.4dBm according to 

(5.3.12) , and the measured power, indicated by the marker in Figure 5.9, 
was about -23.4dBm. The one sided noise spectral density according to 

(5.3.13) was -82dBm/Hz, and the measured value was about -83dBm/Hz 
(-32dBm in the lOOKHz resolution bandwidth). The PLL loop signal to 
noise ratio defined by (5.3.16) was therefore about 30 dB. The average 
cycle slipping time is, according to (5.3.15), over 10 years for SNRl=9.1c 1B 
and B l =380Hz. This implies that the system could achieve PPM word 


Normalized Optimal Threshold Level 
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Figure 5.8. Optimal threshold levels normalized to the peak matched 
filter output vs. detected photons/bit. The parameters used are the 
same as those in Figure 5.7. 
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clock synchronization with virtually no cycle slips at a SNR L much lower 
than 30dB. The system could tolerate at least a 20dB reduction in SNR l 
that might result from higher error rates in the detected PPM pulse 
sequences and/or less frequent occurrences of back-to-back PPM pulse 
pairs. 

To show the effects of misses and false alarms, the threshold level 
applied to the comparator of Figure 2 was changed such that detection 
errors in the PPM pulse sequences consisted mainly misses or false alarms. 
The results of the measurements agreed with those obtained from the 
theoretical analysis, (5.3.18) and (5.3.29). For example, when the threshold 
level was set well above the optimal value such that Pf a «p ms and 
p ms =0.121, the measured signal power and noise spectral density decreased 
by about -2.46 dB and -1.1 dB, respectively. These values were close to the 
values computed from (5.3.18) and (5.3.29), of -2.24dB and -0.90 dB. When 
the threshold level was set well below the optimal value such that 
p ms «Pf a and Pf a =0.0267, the measured noise spectral density increased by 
about 3.0 dB, which again agreed with the value predicted by (5.3.29), of 
3.1 dB. 

Finally, the recovered word clock, the recovered slot clock, and the 
PPM detection/demodulation circuit were all integrated together as a com- 
pletely self synchronized optical PPM receiver. The bit error probabilities 
of the output binary sequences were measured as a function of the average 
number of detected photons per information bit, as shown in Figure 5.10 
(small triangles). The results from the experiments and the numerical com- 
putations under perfect receiver synchronization are also plotted in Figure 
5.10 (small squares and the solid curve). It was found that the performance 
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Figure 5.10. Receiver BER vs. average number of detected signal pho- 
tons per information bit. The small triangles represent the measure- 
ment data when the recovered clocks were used. The small squares 
represent the measurement data with perfect slot and word timing 
recovery. The solid curve are the numerical results assuming perfect 
timing recovery, as given in Chapter 2. The small circles represent the 
measured results when a fixed threshold level of the comparator was 
used. The parameters used were the same as those in Figure 5.7. 
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of the self synchronized receiver was the same as the performance of the 
receiver with perfect synchronization for received optical signal levels 
greater than 25 detected photons/bit. The receiver could acquire and hold 
both slot and word synchronization for received optical signal levels as low 
as 20 detected photons/bit (40 photons per PPM pulse). The small circles 
in Figure 5.10 represent the measured data when the matched filter and the 
comparator were AC coupled and the threshold level of the comparator was 
fixed at zero volts. The measured bit error probabilities were unaffected by 
the use of this close to optimal threshold level except when the received sig- 
nal levels were very low. 
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Chapter 6. Conclusion 

The results of this research showed that a Q=4 PPM direct detection 
optical communication system operating at a source data rate of 25 
Mbits/second using an AlGaAs laser diode transmitter (X=833nm) and a 
silicon APD photodetector could achieve BER<10~ 6 at received optical sig- 
nal energies corresponding to an average of 55 detected photons per infor- 
mation bit. The receiver sensitivity achieved was the closest to the quan- 
tum limit yet reported for an optical communication system of this type. 
Analysis and experiments also showed that reliable PPM slot clock and 
word clock waveforms could be regenerated at the receiver with negligible 
phase jitter at these values of received optical signal energy. The theoreti- 
cal analysis of receiver BER and rms phase errors in the recovered clock 
waveforms agreed very well with the results of actual measurements. 

A new algorithm for computing receiver bit error rates was developed 
which used the nearly exact Webb’s approximation for the APD output 
photocurrent but did not require excessively long computation times. The 
widely used Gaussian approximation for the APD output was found to 
severely underestimate the optimal value of the average APD gain and 
overestimate the receiver BER when the average number of detected back- 
ground noise photons was less than one per PPM slot time. The numerical 
results obtained using Webb’s approximation agreed well with the experi- 
mental measurements under very low background radiation levels. The 
results obtained using the Gaussian approximation converged with those 
obtained with the nearly exact Webb’s approximation when the average 
number of detected photons per slot time due to background radiation noise 
was greater than ten. The mathematical modol nnd the computational 
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procedure developed In Chapter 2 may be applied to other types of optical 
receivers that use APD photodetectors, such as a semiconductor laser radar 
system. 

A transition detector type PPM slot clock recovery scheme was dev- 
ised, analyzed, and implemented. This slot clock recovery scheme was 
demonstrated to be a reliable and simple way to regenerate PPM slot clock 
waveform from the received random PPM signals. A mathematical model 
was developed to accurately compute the rms phase error of the regenerated 
slot clock and it agreed very well with the results of experimental measure- 
ments. The limited bandwidth of the electronics allowed the PPM detec- 
tion circuit to tolerate, without compromising the receiver sensitivity, an 
rms phase error which was less than one percent of a slot clock cycle. The 
experimental slot clock recovery system achieved an rms phase error of well 
below one percent of a slot clock cycle with a PLL loop bandwidth of lKHz, 
or 2X10~ S of the slot clock frequency. Use of the recovered slot clock 
resulted in no observable penalty in the number of detected signal photons 
per bit required to achieve BER over the range of 10" 3 to 10 -7 . The 
theoretical analysis and the experimental measurements also showed that 
the PLL could be replaced with a simple crystal bandpass filter if the 
bandwidth of the filter was one half the PLL noise bandwidth. The only 
disadvantage in using a bandpass filter was the inherent static phase error 
in the output when the input frequency was shifted. 

A novel PPM word timing recovery scheme was devised based on the 
detection of back-to-back PPM pulse pairs in the received random data and 
the use of a PLL. Unlike other existing word synchronization schemes, no 
prescribed synchronization patterns need to be periodically inserted in the 
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data stream and there is no channel use overhead associated with insertion 
of these patterns. Reliable PPM word clock recovery could be obtained as 
long as the transmitted random data were equally likely or back-to-back 
PPM pulses appeared frequently enough. The maximum acquisition time of 
the PPM word synchronization was on the order of a PPM word clock 
cycle. False PPM word synchronization might only occur when the received 
PPM sequence contained too few back-to-back pulse pairs. This could be 
avoided by properly encoding or encrypting the source data. A special PLL 
was devised which synchronized the recovered PPM word clock to the slot 
clock and therefore corrected the random phase error in the recovered word 
clock with respect to the slot clock. The recovered slot clock was more reli- 
able and contained less jitter because the input signal to the slot clock 
recovery system consisted of PPM sequences which had a fixed pulse rate 
regardless of the transmitted data patterns. Use of the recovered PPM word 
clock in the receiver did not cause any measurable degradation in the 
receiver performance when compared to a receiver with perfect word timing. 

The next phase of this work is to increase the binary source data rate 
to 100 megabits/second or more. The assumptions of rectangular pulse 
shape may not be valid because the rise and fall times of the output laser 
light pulses become comparable to the pulse width. The amplifier noise 
spectral density is no longer white but grows with frequency beyond 100 
MHz. Therefore, a more complicated model for the amplifier noise given by 
Personick [32] has to be used. The thermal noise from the APD load resistor 
will increase because the resistance has to be reduced to provide the neces- 
sary preamplifier bandwidth. On the other hand, the PPM slot time over 
which the amplifier noise current is integrated becomes shorter. The overall 



- 130 - 


amplifier noise electron charge integrated over a time slot will be slightly 
larger than at the lower data rates , but the amplifier noise can always be 
overridden by increasing the average APD gain. The major factor which 
limits the receiver performance is the number of detected background noise 
photons per PPM time slot. If the extinction ratio of the laser can be kept 
high, the number of the background noise photons per slot decreases as the 
slot time decreases because the background radiation intensity and the APD 
bulk leakage current are constant. Therefore, the receiver sensitivity is 
expected to improve somewhat at a higher source data rate. 
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Appendix A. Computer Program of BER vs Photons/bit Using 
Webb’s Approximation*. 


INTEGER Q . I ER 

REAL IS, KEEP, nO 
external FF 

common uml , varm.sl , cxI , GNQ , DLO , D0S2 , MMO , NNO , NQO , nO , c 1 3 
c 

c System parameters 

c 

c n1*detectad photons In a slot containing a PPM pulsa 
c nO=datected background nolst photons par slot 
c G=avaraga APO gain 
c 

n1*110 
G*700 . 0 
n0*0.0122 


c 

c 0=numbar of slot In a PPM word 
c R*APO load raslstanca (Ohm) 
c Tr = equ 1 va I ant noise tamparatura (Katvln) 
c Keff=ratlo of Ionisation coefficients of APO 
c Is*APD surface leakage currant In nanoamperes 
c T*PPM slot duration in nanoseconds 
c 

0*4 

R* 1 030 . 0 
Tral 100.0 
KaffaO.OlO 
Is*1 1 .9 
T* 20 . 0 

C 

c Limits and Increments of the summation Indeclas 

c 

c LG* 1 ncramant s of average APO gain 
c 

LG*30 


c 

c MMlalower limit of the outer sum index 
c NN1*uoper limit of the outer sum index 
c Nil*'! ncrements of the outer sum Index 
c 

MMl*n1*l 
rnlareal (nl ) 

fakaff •G*(2.0-t . 0/G ) * ( 1 . 0-kef f ) 
xNN 1 *34 , «f*g*2.»g*rn1 
NN 1 * 1 nt ( xNNI ) 

Nl l*n1 »LG 


c 

c MM0*lower limit of the Inner sum Index 
c NN0*upper limit of the Inner sum index 
c NNO*lncrement of the Inner sum Index 
c 


N00*300 
MMO* 100 
NNO* 1 50000 


* This is the basic version of the Fortran programs used in Chapter 2 for fixed parameter 
values only. Loops may be added to run for several sets of parameter values at a time. 
To use the program, first edit the system parameter section, compile the program, and 
then link it with IMSL (single precision) subroutine. 
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ORIGINAL PAGE is 

OF POOR QUALITY 


e 

c Print out the upper limits »nd Increments to the terminal 
c 

print *. 'N1/0N1/N0/0N0» ' , NN 1 , N 1 1 , NNO , NQO 

C 

c 

c GloPle Constants 
c 

c Constants used In the sums 
c 

0 1 S*FL0AT(n 1 ) *G*G*F 
0 1 = sqrt (D 1 S) 

0 1 S2*2 . *□ 1 S 

OOSanO*G*G*F 

OOssqrt(DOS) 

D0S2*2 . ‘DOS 

0Ll*REAL(n1)*G*P/(F-1.) 

DLO=NO*G*F/(F-1 . ) 

GN1*G«real (nl ) 

GNO=G*nO 
exists *T* 1 .0a- 18 
cx3 = 2 .0*1 . 38e-23»Tr«T*1 .Ce-9/R 
cx4* 1 . 6e- 1 9* I s*T* 1 .e-18 
varm»ex3+ex4 
s1=sqrt(varm) 
c 

c s 1 6*1ntegrat Ion Interval 

c b*sma11ast temr to Pe collected for the outer sum 
c cl3»common factor of part of the Inner sum 
c 

s16*s1*8.0 

b= I . e-9*2 . 0*3 . Ml6*00/real (0-1 ) 
c 1 3*sqrt (2.0*3. 141 6) *00 
c 

c c1*common factor of the sums 
c 

cl 1s>Z.0*3. 14*01 /real (Nil) 
c 1 2 = sart (2.0*3. 14) *00/ (real (q)-1 .0) 
ct*cll»ct2 
c 

c Absolute and relative errors of the integration routine 
c 

AERR*0 . 0 
RERR»0.01 
c 
c 

c Beginning to Calculate 

c 

xx»0 . 0 

do 100 m1»MM1 ,NN1 , N 1 1 

x Israel (Ml )-GN1 

x2»x 1 /OLI+I .0 

x3» 1 . 0/x2 

x5»x1*x1/01S2*x3 

Pmn I »x3*sqrt (x3)*exp(-x5) 

uml areal (ml )*1 . 6e- I 9+Cx 1 

Amt#um1-s IS 

Bml *um1+s 1 6 

c 

c call IMSU numerical Integration subroutine "OCACRE" 
c 

pem1s0CA0RE(FF, Ami ,3ml , A5SR .R ERR , ERROR . IES) 

If (peml.lt.b) goto 150 
100 xx»xx*Pmnl*Peml 
c 

c Print out the actual upper limit of the outer sum index 
c 

150 print ‘.'actual N1» '.ml 


non 
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ORIGINAL. PAGE IS 

OF POOR QUALITY 


c Compute probability of word arror and BER 
c 

PWO« 1 . /C 1 **X 
c2=real ( Q) 
c3=c2/2./(c2-1) 

BER 3 c3»PWD 

c Print out results in file "FOR101.DAT" 
c 

wri t#( 101 , 200) 

200 format( lx . 'MBIT/S N0/T G KEFF IS RTEMP RL0AD ' ) 
c 

c 8R=bit rata 

c Nsphotons/bi t (approximate ns=n1) 

rL*alog(c2)/a1og(2.0) 

L 3 i nt ( rL) 

BR«1 . /T/L* 1 000 . 0 
N»int(N1/L) 

wri teC 101 , 250)BR,n0,G.KEFF f IS.TR.R 
250 format (lx,F5.0,F7.4,F7.0,F5.4,F6.2,2F7.0) 
write( 101 , 270) 

270 f ormat ( 1 x Photons/bi t BER') 

wri te( 101 , 300)N, BER 
300 f ormat (1x.X7.E15. 4) 

900 continue 
stop 
and 
c 

c Subroutine to compute the Integrand 
c 

REAL FUNCTION FF(X) 
real nO 

common uml , varm, s 1 . cx 1 , GNO , DLO , DOS 2 , MMO , NNO, N00 , nO , c 1 3 
c 

c The inner sum from MMO to NNO 
y=*Q . 0 

do 500 mO*MMO . NNO , N00 
umOwrea 1 (MO ) • 1 . 6e- 1 9+cx 1 
y 2 3 (umO-x ) / * 1 

call IMSL subroutine "mdnor" ■ error function 

ca 1 1 mdnor (y2 , y3) 
y4*raal (mO)-GNO 
ySwy4/DLO*1 .0 
y6*i.0/y5 
y 7 3 y4»y4»y6/D0S2 
PmnQwy 6*sqrt CyS)*exp(-y7) 
yy=pmn0*y3 
500 y=yy»y 
c 

c The inner sum from 1 to MMO 
c 

y s*0 . 0 

do 600 m0*1 ,MM0 
umOxreal (mO ) * 1 . 6e- 1 9*cx 1 
y2 3 (um0-x)/sl 
call mdnor(y2,y3) 
y4wreal (mO)-GNO 
y5 3 y4/OL0* 1 . 0 
y 6 s 1 . 0/y5 
y 7=y4»y4«y6/D0S2 
Pmn0 3 y6*sqrt(y6)*exp(-y7) 
yy»pmn0*y3 
600 y»*yy*ys 

y*y»real (n00)+ys 


o n 


-134 


£?'£ WAL PAGE . s 
DF P0 °R QUALITY 


c Th* m«0 tarm 

zl*(cx1-x)/s1 
call mdnor(z1,z2) 
z=z2*c13»«xp(-n0) 
y*y+z 

Tha complata Integrand 
yy I =X-um 1 

yyl-yy l*yyt/2.0/varm 
FF*exp( -yy 2) /a 1 *y 
return 
end 
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